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FOREWORD 


This memorandum presents the results of work performed by Northrop 
Services, Inc. while under contract to the Aero-Astrodynamics Laboratory of 
the Marshall Space Flight Center (NAS8-21810) . This task was conducted in 
response to the requirements of Appendix E-l, Schedule Order No. 3, Technical 
Directive No. 1. Technical Coordination was provided by Mr. Wayne Deaton of 
the Guidance Applications Section (R-AERO-GG) . 
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ABSTRACT 

This memorandum presents a description of the strategy and logic used in 
a space shuttle on-orbit rendezvous targeting program. The program generates 
ascent targeting conditions for boost to insertion into an intermediate parking 
orbit, and generates on-orbit targeting and timeline bases for each maneuver 
to effect rendezvous with a space station. Time of launch is determined so 
as to eliminate any plane change, and all work was performed for a near- 
circular space station orbit. 
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KEY WORDS 


Orbiter - chaser or pursuit vehicle 

Space Station - any target vehicle, satellite 

Shuttle Launch Vehicle - booster plus orbiter configuration 

Intermediate Orbit - a phasing orbit for the orbiter on-orbit used to alleviate 
large phasing differences between vehicles. (- 100 n mi Parking orbit for this 
analysis) 

Constant Delta Height (CDH) - a height differential existing between the 
orbiter and the space station (an orbit approximately 10 n mi below or above 
the space station). Same as coelliptic orbit. 

Transfer Phase Initiation (TPI) - A point on the CDH orbit when gross rendez- 
vous conditions have been met in order to make the final transfer to the 
rendezvous point. 

On-Orbit - the pursuit vehicle after insertion and before rendezvous during 
all of its intermediate phasing orbits. 
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SYMBOLS 

Wedge angle between planes at TPI 

Range angle difference at TPI after isolation 

Nodal difference at TPI 

Lift-off time correction to compensate for nodal regression 

First pass Range angle difference using first guess two-body 
targeting 

Space fixed launch coordinate system 
Launch aximuth 

Geodetic latitude of launch site (28.608°) 

Sun vector right ascension 
Sun vector declination 

Universal Time -measured from midnight Greenwich to launch meridian 
longitude of launch site 

difference in range angles of orbiter and space station at time 
of orbiter insertion 

insertion latitude, a function of inclination of the space station 
desired inclination for targeting purposes 

Semi-latus rectum of CDH orbit 

eccentricity of CDH orbit 

earths rotational velocity 

Semi-latus rectum of orbiter on-orbit during Hohmann transfer 

eccentricity of orbiter on-orbit during Hohmann transfer 

Universal time of lift-off 
Range angle 

true anomaly 
argument of perigee 
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SYMBOLS (Continued) 


DS 


DN 


LS 

x t ,x t 

T, 


Semi-major axis 
eccentricity 

desired insertion latitude for southerly launch 
desired insertion latitude for northerly launch 
Range angle at the desired latitude 
State vector position and velocity of orbiter 
State vector position and velocity of space station 
Time of orbit insertion 
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instantaneous angle from the suns projection vector on orbital 
plane to the TP I point 

same as above but is the desired input value 

Unit vector in -the equatorial plane and through the launch 
longitude 

descending node of space station referenced from space-fixed 
shuttle launch meridian in the equatorial plane 

descending node of the orbiter referenced from space-fixed shuttle 
launch meridian in the equatorial plane 

desired range angle difference between vehicles at TPI 

difference between actual and desired range angle difference, this 
value to be driven < .05 in the isolation logic 

Range angle of space station measured from the descending node 
w.r.t. equatorial plane 

Range angle of orbiter measured from descending node w.r.t. 
equatorial plane 

Range angle of the space station measured from the common 
(ascending) node of the space station and the orbiter planes 

Range angle of the orbiter measured from the common (ascending) 
node of the space station and the orbiter planes 

difference in the range angles of the space station and the 
orbiter 
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WATP 

WATOL 

TSTI 


SYMBOLS (Concluded) 

Wedge angle between the space station and the orbiters plane 

Tolerance to select which &<j> to use (for example, WATOL < .1 : : 
A+ = + T ~ Acf> = 4 > nt - <fr OT ) 

Time of Circularization 
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Section I 

INTRODUCTION 

This memorandum is primarily an equation defining document containing 
the basic targeting equations in flowchart form to create targeting conditions 
at lift-off for the shuttle launch vehicle. This also includes the method of 
determining the on-orbit timeline of thursting events* during orbital maneu- 
vers and also determines the Universal Time of lift-off. 

The basic mission profile considered for this targeting procedure includes 
boost to insertion and three impulsive maneuvers, as listed below, to establish 
a constant delta height position (Figure 1-1). 

• Insertion (50 x 100 n mi) 

• Circularization at apogee (100 n mi) 

• Perigee impulse (- 100 x 265 n mi) 

• Coelliptic impulse (z 260 n mi) 

The launch azimuth (A^) , inclination (i) and node (0 N ) for the launch 
phase are generated to achieve orbiter/satellite rendezvous. These are 
generated in such a manner as to achieve orbiter/satellite rendezvous with 
coplanar conditions near rendezvous and with the proper phase and coelliptic 
height differential at TPI. 

The given task assignment was to build a space shuttle on-orbit rendez- 
vous targeting computer program that would depend only upon a target satellite 
ephemeris and the initial in-plane orbital conditions of the space shuttle 
(50 x 100 n mi). The computer program was to establish lift-off time for the 
space shuttle so as to require no plane change in the ascent portion of flight, 
or on-orbit portion of the rendezvous mission. The computer program estab- 
lishes a timeline of the thrusting events and guidance targeting requirements. 

*This targeting procedure is developed with impulsive maneuver simulations. 
Using these targeting values on-orbit will result in ignition time deviations 
for each maneuver. This could be alleviated by simulating finite bums with 
the targeting deck itself. 
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Figure 1-1. COPLANAR PROFILE DEPICTING TIME BASES 
FOR NEAR-CIRCULAR RENDEZVOUS 
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Care was to be taken to minimize the number of instructions and storage 
requirements of the program so that it would be possible to have an on-board 
shuttle rendezvous capability. The Coordinators flowcharts were to be used, 
and deviations were to be made whenever necessary and storage instruction 
could be reduced. 
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Section II 
DISCUSSION 

The shuttle, being a performance critical vehicle, should be targeted to 
a zero plane change, on-time ascent to orbit (50 x 100 n mi) flight profile 
(as well as to basic satellite delivery missions) . The shuttle should not be 
burdened with a requirement for a rendezvous launch window since this would 
degrade the payload delivery capability. The procedure presented here will 
allow launches to be achieved at each in-plane point. One in-plane point will 
occur for a northerly launch opportunity and the other for a southerly launch 
opportunity. These conditions occur twice per day, 365 days/year. These two 
launch opportunities that occur each day are only restricted if the launch 
site is too close to the in-plane point to allow pre-flight analysis to be 
performed before the launch. With more restrictive launch vehicles (short 
systems liftimes) , the correct in-phase and in-plane condition (rendezvous 
compatible) has to exist to achieve a rendezvous; but, this is not a require- 
ment for the targeting technique presented in this memorandum. An intermediate 
near-circular phasing orbit at the apogee of the shuttle launch vehicle 
50 x 100 n mi insertion orbit will eliminate the space station in-phase 
requirement at orbital insertion. (If the relative catch-up rate between the 
100 n mi intermediate phasing orbit and the space station is not sufficient 
to null out phase differences, the use of an intermediate stay orbit at a 
higher altitude will be necessary.) An intermediate phasing orbit exists so 
that phase angle differences between the two vehicles can be eliminated by 
exploiting the difference in their respective orbital periods. 

Other advantages of this targeting technique include: 

• Launch vehicle performance variations will merely change the range 
correction of the terminal rendezvous maneuvers without causing 
unacceptable performance losses. 

• Eliminates high closing rates of the orbiter w.r.t. the space 
station, which might be encountered when using direct rendezvous 
techniques and their resulting performance losses. 
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@ This technique allows launch opportunities to occur on a daily basis 

without degrading the payload delivery capabilities. This is important, 
for example, when considering the shuttle launch vehicle configuration 
which requires many launches each year for economical reasons. 

& If count-down is delayed the next opportunity can be utilized. 

The targeting program generates complete targeting based upon space 
station ephemeris data. This is accomplished by assuming that the Manned 
Space Flight Tracking Network (MSFN) has made available the epoch (Universal) 
time when the launch site will be contained in the space station plane, based 
upon spherical trigonometry and also the ephemeris at this time. The orbital 
elements (node, inclination, eccentricty, etc.) describing the position of 
the space station at the in-plane time (U.T.) are presented in Section V. 

The periodic perturbations of the space stations' inclination were 
determined and accounted for in the targeting procedure by using a rapid 
integration algorithm to advance the space station to the insertion latitude 
of the shuttle (at present a variable step size Runge Kutta numerical 
integration scheme is utilized) . 

The effects of orbital nodal regression are corrected by adjusting the 
shuttle launch vehicle lift-off time while maintaining the same ascent 
targeting parameters. The amount of nodal regression depends on the transfer 
orbits necessary to satisfy phasing requirements, navigation update require- 
ments and lighting requirements. 

The ascent trajectory was programmed as a functional representation 
of an ascent profile. This is presently a sixth order curve fit polynomial 
as shown in the flowchart on page D-6. Future work in this area includes 
curve fit techniques using exponential curves and other types of fits which 
will improve curve-fit accuracy and reduce empherical curve-fit coefficients. 
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Section III 

RENDEZVOUS TARGETING TECHNIQUES 

The procedures for effecting rendezvous includes integration of the 
space station to the insertion latitude (ijO to determine the desired inclina- 
tion (I ) for the orbiter insertion. This causes the orbiter to have the 
same mean inclination as the space station at insertion. This procedure is 
necessary to account for periodic variations in the inclination of the 
space station orbit about the oblate earth. The variation of inclination 
versus time from insertion and time after circularization for both vehicles 
is presented in Figures 3-1 and 3-2. These figures depict variations with 
approximately the same mean inclination. Similar results, at a point on- 
orbit after the apsidal rotation maneuver where the orbiter is phasing 10 n mi 
below the space station, are presented in Figure 3-3. As can be observed at 
this point, the variations in inclinations are almost in-phase and thus nearly 
synchronized. This is desirable for rendezvous targeting to alleviate unnec- 
essary plane change during coelliptic coast. 

The desired inclination (1^) for targeting purposes dictates the ascent 
targeting parameters for the shuttle booster /orbiter launch configuration. As 
shown on page D-6 of the flowchart the insertion conditions for the southerly 
and northerly launches are a function of the desired inclination. The launch 
azimuth, descending node, insertion time, and range angle are presently least 
square curve fit functions of the desired inclination (L) . 

A quick-look two-body analysis of the on-orbit phasing is executed after 
orbit insertion. Many of the two-body parameters (page D-8) are used for the 
initialization of the isolation technique for its "first guess". 

The time bases for each on-orbit maneuver are given in Table 3-1. These 
time bases occur approximately 200 seconds prior to the actual maneuver. The 
actual times will be presented in Section IV. This timeline includes the 
insertion time, the time of circularization, time of perigee burn out of the 
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Figure 3-1. ORBITER AFTER INSERTION WITH APPROXIMATE SAME MEAN INCLINATION 
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Figure 3-2. ORBITER AFTER CIRCULARIZATION WITH APPROXIMATE SAME MEAN INCLINATION 
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Figure 3-3. SHUTTLE RENDEZVOUS INPLANE LAUNCH GEOMETRY 
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100 n mi intermediate orbit, time of coelliptic maneuver, and the time of 
Transfer Phase Initiation (TPI). 

Table 3-1. TIME BASES FOR PREPARATION OF MANEUVERS 


SYMBOL 

DEFINITION 

UNIT 

T 1 

Insertion time of orbiter 

sec 

TTEST 1 

Preparation for the circularization 
maneuver begins at this time 
(approximately 200 seconds before apogee 
of initial insertion orbit of the 
orbiter) 

sec 

TTEST 2 

Preparation for perigee maneuver out of 
the near-circular phasing orbit begins 
at this time 

sec 

TTEST 3 

Preparation for the coelliptic maneuver 
begins at this time 

sec 

TSBST 

Time for the transfer phase initiation 
(TPI) 

sec 


The mission time of insertion and time of circularization will not be 
changed during the isolation loop. This is a constraint which must be met 
to maintain the same ascent targeting parameters. 

The time of perigee burn out of the 100 n mi intermediate orbit is a 
variable and depends on the desired value of the phase angle difference (A<j> D ) 
at TPI. It also depends on the oblate earth effects on the "first guess" 
two-body timing. 

The effects of orbital nodal regression are corrected by adjusting the 
shuttle lift-off time. Once the desired phase angle ( Acp at TPI is isolated, 
then the existing nodal error is mapped into a Universal Time correction at 
lift-off. This is demonstrated on page D-18 of the flowchart. 
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The manner in which the sun's declination and right ascension are 
evaluated is illustrated in Figure 3-4, and the suns position in the launch 
coordinate system is depicted in Figure 3-5. Knowledge of the suns position 
is necessary in the targeting procedures when proper lighting is considered. 

A general flowchart of the rendezvous targeting technique is presented in 
Figure 3-6. Detailed flow of this targeting procedure is included on pages 
D-2 through D-20. 

A typical mission profile is illustrated in Figure 1-1. The orbiter is 
inserted into a 50 x 100 n mi orbit. A coast to apogee occurs where an 
apogee burn is made to circularize into the 100 n mi circular orbit (TEST1) . 
After circularization a coast of at least a half-orbit is necessary (and is 
handled by the scale factor input SFN01) . A value of SFN01=0.5 insures at 
least a half-orbit before the perigee burn onto a Hohmann transfer at time 
TTEST2 . This scale factor can be initialized to any desired value. More 
stay time would be desired if phasing or lighting constraint is to be satis- 
fied. The purpose of extra stay time would be to insure time needed for real 
time preparations. The scale factor for the Hohmann phasing (SFN02) and the 
coelliptic phasing, 10 n mi below or above target, (SFN03) will insure extra 
stay time in all phasing orbits until rendezvous is accomplished. This extra 
stay time will enforce adequate time for crew and orbiter check-out, orbit 
evaluation and system checkout, propulsion checkout, tracking acquisition, and 
navigation up-date. Any realistic targeting technique has to provide this 
extra controlled stay time for real time targeting. 

After the perigee maneuver at time base TTEST2, a coast of approximately 
a half an orbit brings the vehicle to an intersection with the coelliptic 
orbit. The derivation of the equations for determining the intersection of 
the near-Hohmann transfer with the Constant Delta Height (CDH) orbit at time 
base TTEST3 is presented in Appendix A. The equations necessary for deter- 
mining the desired values for the differential height are included in Appendix 
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Figure 3-4. RIGHT ASCENTION AND DECLINATION OF SUN 
WITH RESPECT TO LAUNCH MERIDIAN 
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Figure 3-5. ROTATIONS FROM LAUNCH COORDINATE SYSTEM TO SUN VECTOR 
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SHUTTLE ASCENT TARGETING (Concluded) 
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B. These equations have been presented in a earlier publication (ref. 1), 
but in a different manner. 

A pictorial illustration of the position of the orbiter and space station 
at insertion (T^ + T^) and at lift-off (T^) for the tracking network is given 
in Figure 3-7. The A<j> angle represents the range angle difference between the 
space station and the orbiter at insertion time T^ . It is reasoned that the 
tracking network will supply the ephemeris of the target when it is in-plane 
at U.T. of T n or Tg, and not the ephemeris that the target vehicle has at any 
acquisition time T. The ephemeris will be used to determine the time devia- 
tion (AT^) for lift-off from this in-plane time (T^) which will result in 
coplanar on-orbit phasing near the rendezvous point. It is not a necessary 
criterion, as stated earlier, that this be a rendezvous compatible orbit, so 
any A<j> relation may exist at lift-of f/insertion and rendezvous can be 
accomplished through proper on-orbit phasing. 


1. Wesselj V. W. j Bentley , E. L. 3 and Sport 3 R. H. 3 "Guidance Equations for 
AAP-4 S-IVB/LM-ATM Unmanned Rendezvous" 3 Northr op-Hunt sville Technical 
Report TR-795-9-531 , April 1969. 
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Section IV 

RESULTS AND CONCLUSIONS 

The rendezvous targeting program was developed to generate targeting 
conditions for the shuttle launch vehicle at launch. The desired inclination 
(I ) and launch azimuth (A ) at lift-off can be determined to achieve rendez- 

U Lt 

vous with near-circular target satellites at various inclinations and various 
altitudes. Also, the time of launch (Universal Time, U.T.) and the timeline 
bases from lift-off have been determined for the orbital maneuver to accomplish 
rendezvous. 

Verification of the targeting scheme included a total of 30 cases being 
run with varying phase relationships of the space station (0 < A<j> < 2ir) at the 
time of orbiter insertion (Figure 4-1). Included were cases with lighting 
constraints, northerly and southerly launch opportunities, and different phase 
relationships at transfer phase initiation. 

Several cases were run for a northerly launch, without a lighting con- 
straint. The phase relation of the orbiter at. TPI is below and behind the 
space station by a 10 n mi height differential (the orbiter lags the space 
station by a desired A 4> = -0.29 degree). Isolated phase angles at TPI of 0.29, 
0.32, 0.28, and 0.24 degree were obtained (as shown in Table 4-1), which are 
all within the desired tolerance of 0.05 degree. Also, the inplane conditions 
at TPI are within acceptable limits as can be observed from the values of Ai 
and A4> n . These inplane conditions could be improved, if desired, by decreasing 
the tolerance (of 0.02 degree) on pages D-18 through D-20 of the flowchart. 
Also, the timeline for the on-orbit maneuvers is listed in Table 4-1 for dif- 
ferent range angles of the space station at the time of orbiter insertion 
(255, 345, 75, and 165 degrees). The adjustment required in the lift-off 
time is listed as AT L q, with the negative values representing launch before 
the spherical in-plane point. 

Similar results for a southerly launch with a lighting constraint are 
presented in Table 4-2. The desired sun angle input was 110.0 degrees. 
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Figure 4-1. INCLINATION SYNCHRONIZATION DURING COELLIPTIC COAST -ON-ORBIT DECK- 
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Table 4-1. LIFT-OFF TIME CORRECTION FOR NEAR-CIRCULAR TARGET ORBITS 
FOR NORTHERLY LAUNCH WITH NO LIGHTING CONSTRAINT 


V '\ V JTRUE ANOMALY 

45° 

135° 

255° 

315° 

ABOVE ^\ R ANGE 
AND \.ANGLE 

AHEAD 

255° 

345° 

75° 

165° 

Insertion 

0 hr. 6 m. 
11 s. 

0 hr. 6 m. 
11 s. 

0 hr.. 6 m. 
11s. 

0 hr. 6m 
1.1 s. 

Circularization 

0 hr. 48 m. 
27 s. 

0 hr. 48 m. 
27 s. 

0 hr. 48 m. 
27 s. 

0 hr. 48 m. 
27 s. 

Peri gee 

1 hr. 26 m. 
7 s. 

7 hr. 14 m. 
41 s. 

12 hr. 23 m. 
5 s. 

18 hr. 36 m. 
11 s. 

Constant Delta Height 

2 hr. 14 m. 
33 s. 

8 hr. 4 m. 
42 s. 

13 hr. 13 n. 
17 s. 

19 hrs. 27 m. 
3 s. 

Transfer Phase Initiation 

5 hr. 48 m. 
19 s. 

11 hr. 10 m. 
59 s. 

16 hr. 37 m. 
46 s. 

21 hr. 58 m. 
17 s. 

Sun Angle (deg) 

N. A. 

N. A. 

N. A. 

N. A. 

Ai (deg) 

5.3 x 10' 4 

1.1 x 10' 2 

6.9 x 10' 4 

4.4 x 10' 4 

A4> (deg) 

0.29 

0.32 

0.28 

0.24 

aQ n (deg) 

-4.0 x 10' 4 

-1 .3 x 10' 2 

-2.8 x 10" 4 

4.6 x 10‘ 4 

aT lq (sec) 

-284.7 

-226.4 

-199.06 

-143.6 

ai> tb (deg) 

-.81 

-2.21 

-4.43 

5.94 

SFN03 (unitless) 

1.5 

1.5 . 

1 .5 

1.5 
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Table 4-2. LIFT-OFF TIME CORRECTION FOR NEAR-CIRCULAR TARGET ORBITS 
FOR SOUTHERLY LAUNCH WITH LIGHTING CONSTRAINT 


TRUE ANOMALY 



RANGE 
\ ANGLE 


Insertion 

Ci rculari zation 

Perigee 


Sun Angle (deg) 
Ai (deg) 
ao (deg) 
a9 n (deg) 
aT lq (sec) 
ao tb (deg) 

SFN03 (unitless) 


45° 

135° 

225° 

315° 

255° 

345° 

75° 

165° 

0 hr. 6 m. 
11 s. 

0 hr. 6 m. 
1111 s. 

0 hr. 6 m. 
11 s. 

0 hr. 6 m 
11 s. 

0 hr. 50 m. 
48 s. 

0 hr. 60 m. 
4 48 s. 

0 hr. 50 m. 
. 48 s. 

0 hr. 50 m. 
48 s. 

18 hr. 27 m. 
36 s. 

2 hr. 39 m. 
29 s. 

7 hr. 51 m. 
52 s. 

13 hr. 30 m. 
17 s. 

19 hr. 16 m. 
59 s. 

3 hr. 30 m. 
3 s. 

8 hr. 42 m. 
23 s. 

14 hr. 19 m. 
59 s. 

26 hr. 22 m. 
50 s. 

8 hr. 37 m. 
27 s. 

14 hr. 33 m. 
22 s. 

18 hr. 50 m. 
34 s . 

109.91 

109.88 

109.92 

109.84 

8.17 x 10' 4 

2.8 x 10' 3 

7.9 x 10" 3 

1 .27 x 10" 3 

-.332 

-.256 

-.317 

-.278 

-1.87 x 10' 4 

-2.53 x 10‘ 3 

-9.3 x 10' 3 

5.2 x 10' 4 

319.65 

190.1 

235.9 

275.42 

35.5 

13.36 

22.7 

16.11 

4.0 

3.0 

3.5 

2.5 
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The results presented in Tables 4-1 and 4-2 are for effecting rendezvous 
with the baseline target in a 270 n mi orbit with an approximate 55 degree 
inclination. Rendezvous targeting was accomplished in all cases considered, 
with a time constraint of approximately 24 hours to the TPI point. 

It should be noted that rendezvous with satellites at altitudes other 
than 270 will result in violation of a 24 hour time constraint, but targeting 
is still possible. This violation will most likely happen when the target 
satellite has a lower altitude and thus additional phasing in the 100 n mi 
phasing orbit will be required to alleviate large phase differences which may 
exist. 

The executed listing presented as an example in Appendix C gives the 

eccentricity vector e, angular momentum (h) and delta velocity required 

(AV ) at each maneuver time to effect each burn. These on-orbit targeting 
R 

conditions at each maneuver time can be used as inputs for any guidance 
package to simulate that particular orbital maneuver. 

Results to date show that the present rendezvous targeting deck will 
establish lift-off time and on-orbit targeting parameters to effect gross 
rendezvous at TPI. 
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Section V 

PROGRAM INPUTS AND OUTPUTS 

5.1 INPUT 

The rendezvous targeting deck was programmed in Fortran IV language for 
use on the CDC-3200 computer. Inputs to the program are described in the 
following text, and are listed in Table 5-1. 


Table 5-1. INPUTS FOR NEAR-CIRCULAR RENDEZVOUS 


n 

FORTRAN 

ALFA-NUMERIC 

NAME 

DEFINITION OF SYMBOL; 

SCALE FACTOR, FLAG, VARIABLE, ETC. 

UNITS 

ISURF 

ISURF 

=1 , Boost cut-off surface 
=0, Steady state trajectory comp, 
(see page D-6 of flowchart) 

Uni tl ess 

ILIG 

ILIG 

=1, Lighting constraint considered 
=0, No lighting considered 

(see page D-10 of flowchart) 

Unitless 

V A 6 

A ( 7 ) 

Polynomial coefficients as a function 
of inclination to determine latitude 
of insertion for northerly launch. 

Deg 

B (T B 6 

B(7) 

(Same as above for southerly launch). 

Deg 

V C 6 

C(7) 

Polynomial coefficients as a function 
of inclination to determine azimuth 
of insertion for northerly launch. 

Deg 

V°6 

0(7) 

(Same as above for northerly node). 

Deg 

E (f E 6 

E(7) 

(Same as above for northerly 
time-of-insertion) . 

Deg 

F cT F 6 

F (7 ) 

Polynomial coefficients as a function 
of inclination to determine azimuth 
of insertion for southerly launch. 

Deg 

V G 6 

6(7) 

(Same as above for southerly node) 

Deg 

V H 6 

H (7 ) 

(Same as above for southerly 
time-of-insertion). 

Deg 

V Q 6 

Q(7) 

Range angle of insertion (northerly). 

Deg 

V S 6 

S ( 7 ) 

Range angle of insertion (southerly). 

Deg 

T 

T 

The universal time ephemeris data 

received from tracking station. 

• 

Sec 
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Table 5-1. INPUTS FOR NEAR-CIRCULAR RENDEZVOUS (Continued) 


MATH 

SYMBOL 

FORTRAN 

ALFA-NUMERIC 

NAME 

DEFINITION OF SYMBOL; 

SCALE FACTOR, FLAG, VARIABLE, ETC. 

UNITS 

t n 

TN 

The universal time of the in-plane 
opportunity for northerly launch. 

Sec 

T S 

TS 

The universal time of the in-plane 
opportunity for southerly launch 

Sec 

TTOL 

TTOL 

Maximum time necessary to perform 
pre-flight analysis using this 
targeting deck. 

Sec 

HAP 

HAP 

Altitude of apogee of orbiter 
insertion ellipse 

N MI 

h per 

HPER 

Altitude of perigee of orbiter 
insertion. 

N MI. 

a n 

AN 

Semi -major axis of space station 
received from tracking network for 
northerly opportunity (T^). 

M 

e N 

EN 

Eccentricity of space station received 
from tracking network for northerly 
opportunity (T^). 

Unitless 

1 N 

X ENCN 

Inclination of space station received 
from tracking network for northerly 
opportunity (T^). 

Deg 

0 NN 

TNNN 

Descending node for northerly launch 
(TN). 

Deg 

“PLN 

ALFAN 

Argument of perigee for northerly 
launch (measured from descending node 
opposite direction of flight). 

Deg 

*N 

PNIN 

True anomally of space station for 
northerly launch. 

Deg 

A S 

AS 

Semi -major axis for southerly launch. 

M 

e S 

ES 

Eccentricity for southerly launch 

Unitless 

n s 

XENCS 

Inclination for southerly launch. 

Deg 

0 NS 

THNS 

Node for southerly launch. 

Deg 

a PLS 

ALFAS 

Argument of perigee for southerly 
launch . 

Deg 

*S 

PHIS 

True anomally for southerly launch. 

Deg 

*L 

PHI 

Geodetic latitude of launch site 
measured from equatorial plane. 

Deg 
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Table 5-1. INPUTS FOR NEAR-CIRCULAR RENDEZVOUS (Concluded) 


MATH 

SYMBOL 

FORTRAN 

ALFA-NUMERIC 

NAME 

DEFINITION OF SYMBOL; 

SCALE FACTOR, FLAG, VARIABLE, ETC. 


»L 

XLAMAL 

Longitude of launch site measured 
negative west of prime meridian. 

Deg 

6 svd 

BSVD 

Desired sun angle, measured from 
sun projection on space-station 
plane in direction of flight to TPI. 

Deg 

a, b, c 

Al, Bl, Cl 

Coefficients for calculation of the 
declination angle of the sun W.R.T. 
in the equatorial plane, (see page 
D-10 of the flowchart) 

Unitless 

e T0L 

TOLE 

When the space station gets within 
this tolerance, a simplified logic 
for the space station in a circular 
orbit will be inacted (eccentricity 
tolerance) . 

Unitless 

t y 

ty 

Number of days past January 1 of 
launch year. 

Days 

aH d 

DLHD 

The desired differential height for 
the orbiter: > 0 :: CHD below target; 
<0 :: CDH above target. 

N MI 

AHg 

DLHB 

A bias used to insure that the trans- 
fer orbit will intersect the C.D.H. 
orbit. 

N MI 

SFN01 

SFN01 

Scale factor for the initial orbiter 
insertion orbit. (Generally = .5) 

Unitless 

SFN02 

SFN02 

Transfer orbit scale factor for inter- 
mediate phasing orbit (SFN02 = .5 for 
second orbital intersection). 

Unitless 

SFN03 

SFN03 

Scale factor for phasing time in the 
coelliptic C.D.H orbit (normally = 
1.5). 

Unitless 

SLM 

SLM 

Slope of the A<j> = f(AH) curve assumed 
to be linear. 

Unitless 


The first input card contains two fixed point options with a 212 format. 
Presently the first option ISURF is flaged as 1. This designates that a 
sixth order polynomial curve fit will be utilized for describing the Shuttles 
insertion surface (Figure 3-3). A future mode may be programmed to execute 
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a steady state trajectory. When this mode is developed the user would read 
ISURF=0. The second option ILIG is for the lighting constraint. If ILIG=1, 
lighting is considered and future inputs will include a * c » aS 

described in the input nomenclature. 

The format for the remaining inputs is 6E13.8. Aq through Sg contain 
the coefficients for the curve fit surface of the orbiters cut-off. These 
are contained on the next 20 cards . 

Input on card 22 are the universal times from the tracking station, 
along with the radius of apogee and perigee of the orbiter insertion orbit. 
Card 23 provides input for the ephemeris for the space station at the time 
(U.T.) the launch site is in-plane with the space station for a northerly 
launch opportunity. Similar values for the southerly launch opportunity are 
input on card 24. The latitude of the launch site longitude desired 
sun angle Bg VD « and coefficients for calculation of the suns declination 
Al, Bl, Cl are input on card 25. Cards 26 and 27 will be changed by the user 
as different mission profiles are desired. These cards contain the desired 
differential height (AH) for the final phasing orbit (coelliptic) before TPI. 
The desired phase angle (A<j>) at TPI is determined as a function of AH and is 
presently read in as a linear function with a slope SLM. 

Three flags are input which represent whole or fractional stay time 
periods in each of the orbiters phasing orbits. SFN01 and SFN02 will be input 
and will stay fixed. SFN03 can and will be "bumped" if the isolation results 
in orbit coast periods in the coelliptic orbit is less than SFN03 times the 
orbital period. That is, when the stay time in final coelliptic orbit between 
the constant delta height maneuver and the TPI maneuver is less than SFN03 
orbits (Note Page D-14 of the flowcharts), then SFN03 will be bumped by .5 
and reinitialized. 


A list of sample input data is presented in Table 5-2. It should be 
noted that only two coefficients are listed for each surface or polynomial 
curve fit variable. 
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5.2 OUTPUT 

The sample output (executed listing) presented in Appendix C is for a 
northerly launch opportunity with lighting considered. The first two pages 
yield the input values from the tracking station and the "first guess" two- 
body analysis of the total mission. The last variable printed on the second 
page, DVIT, gives the delta velocity budget requirement for all on-orbit 
maneuvers; but, does not include the values for TPI and TPF. 

The listing has comment cards throughout, describing each maneuver, and 
gives on-orbit targeting requirements (e,h,Av). Both Universal Time and 
mission time from lift-off for each maneuver is located at the top of the 
page, along with the state variables of the orbiter and space station. 

The last two pages present the final isolated values at the TPI point; for 
example, sun angle, A<j>, Ai, A0 N> and, also, the state vector of the space 
station in the updated coordinate system at the time of lift-off and orbit 
insertion. The very last print statement yields the updated time-of -launch. 
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Appendix A 

INTERSECTION OF NEAR-HOHMANN TRANSFER WITH CDH ORBIT 

A maneuver at the second orbital intersection of the transfer eclipse 
with the CDH orbit will place the orbiter coelliptic with the space station. 
Thus, a method had to be determined to compute the true anomaly of the orbiter 
at the desired second orbital intersection. A solution to this problem is 
possible if the two-body polar equations for position of each orbit are 
equated and then solved for the true anomaly of the intersection. The deriva- 
tion for determining the intersection point follows. 

Considering the equation 

Aa = a T - a p 

where a T is the argument of perigee of the space station orbit and a p is the 
argument of perigee of the orbiter orbit, then 

0 S = 0 p + Aa 

where 0 is the true anomaly of the CDH orbit and 0 is the true anomaly of 

J * 

the orbiter at the intersection point. 

Then, equating the position equations. 



1 + e p cos 0 p 1 + e g cos (0 p +Aa) 


or, 

P P + e S P P cos ^ 9 p + Aa ) = P s + e P p s cos 9 P 

and 

e g P p cos (0 p + Aa) - e p P g cos 0 p = P g - P p 
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Making use of the trigometric identity of the cosine of the Siam of two angles, 
e g P p {cos 6 p cos Aa - sin 0 p sin Aa} - e p P g cos 0 p = P g - P p 
Factoring out cos 0p: 

sin 0 p (-e g p p Aa) + cos 0 p (e g P p cos Aa - e p P g ) = P g - P p 

Now let 

8 = -e g P p sin Aa 
A = e g P p cos Aa - e p P g 

P - p _ p 

o S P 

then; 

6 sin 0„ + A cos 0„ = P 
P P o 

A cos 0_, = P - 8 sin ©„ 

P o P 

A 2 cos 2 0 = P 2 - 2P 8 sin 0„ + 8 2 sin 2 ©„ 

Poo P P 

2 2 
cos 0 p = 1 - sin 0 p 

A 2 (l - sin 2 0 ) = P 2 - 2P 8 sin 8_ + 8^ sin 2 0 
Poo P P 

2 2 2 2 2 2 
A - A sin 0 = P - 2P 8 sin 0 D + 8 sin 0„ 

P o o P P 

and 

(-8^-A ) sin 2 0 + 2P 8 sin 0 + A 2 - P =0 

P o P o 
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N-1 0&2 


In order to sove this quadratic, let 


2 2 
a = -e - a 


B = 2P 6 
o 


2 2 
C = A - P 

o 


and the equation is solved by 

-B + A 2 ~ 4AC 


sin 0 p = 


2A 


This equation is derived as a sine function instead of a cosine function 
as in reference 1. The sine function is positive in the second quadrant and 
negative in the third quadrant. The solution that is negative should be 
selected and placed in the third quadrant (since the transfer is a near- 
Hohmann). This will always select the second orbital intersection, that is, 
select sin 0 p < 0 : : 0 p = -ir - sin (0 p ) 


A-3 



NORTHROP SERVICES. INC. 


M-1082 


Appendix B 

CONSTANT DELTA HEIGHT IMPULSE 

The delta velocity for the impulse into the CDH orbit below or above the 
space station is computed using two-body equations. Forcing the CDH orbit to 
be coelliptical with the space station can only be achieved by having the same 
differential height (AH) at apogee and perigee. Thus, to insure the AH will 
be the same at apogee and perigee, the following equation was developed (see 
reference 1 for complete derivation): 

2 RRP 

AH + (RRP-RAT-RPT ) AH + RPT • RAT + — (RPT-RAT) 

RRP 

cos e D - (RAT+RPT) = 0 


Letting 

A = 1 

B = RRP-RAT-RPT 


C = RAT • RPT + 


RRP 


cos 0 D • (RPT-RAT) - 


RRP 


(RAT+RPT) 


Then 

A(AH) 2 + B (AH) + C = 0 

and 

, TT _ -B ± A 2 - 4AC 


If the coelliptic CDH orbit is above the space station then B = -B. 


These equations are incorporated into the logic on page D-14 of the 
flowchart as can be observed from this flowchart, once AH is computed it is 
utilized to construct the conic parameters of the CDH orbit. 
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Appendix C 

SAMPLE OUTPUT: NORTHERLY LAUNCH 


c-i 



TIME J.OSOBOOnO' ; jJ AT 6.3792069CE 06 ET 1 . 000000009-05 XExCTO 5.50000000E 01 

ThMO 1.57 ji:i: 0 i 5: :-2 ALEato 1 .50000000E 02 P M I 1 0 1.350000005 02 
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PARAMETERS rlJP 50/1J0 N.M. PHASING 0R8IT/ 


RHP 6.47(]76600C 06 
AP 6.51706600C 06 


VVP 7 . 8764368.E 03 GAMHAP-2 , 81 0458 J4E -1 0 EP 7 . 10442689E-0J 

Hap 1.000000UE 6? HPP 4.999999026 01 PMDOTP 6 . B7J66 714E-02 


parameters for the target orbit 


RRT 6.87755324= 06 
AT 6.87879782E 06 


VVT 7.61364626E 03 GAmATO-2 . B5674250E-02 ET 5 . 30393840E-04 
HAT 2.72289576E 02 HPT 2.68349542E 02 PHOOTT 6 . 34051019E-02 


CATC" uP rate a NO ANGLE F hr the half OROIT OF the 50X100 nn phasing orbit 


ORBITAL CATCH UP RATfc. 5. J5156949F.03 
ANGLE OF CATCH UP r 1 .'401 U0222P (!1 

TIHE AT APOr.EE OF 50X100 ORBIT* 2.9698H35E 03 
FIRST maUuVFR TO CIRCULARIZE «o<100 AT ITS APOGEE/ 

RLP 6.567366020 u 6 VCP 7.79304310E 33 TAUCP 5.29175108E 03 PDOTCO 6 . B0304097E-02 
DL°0 2 8. n 7 ?6 8494"-o6 ULMR ? 0 T-22379887E 31 DELV2 2-7731B95JE oi T 3 5-63568689E o3 


Second. burn tra! sfer out of 100 nm circular towards coflliptic 


orbital HE R I Ops 5..74601 73E 03 
mean OkBITai rate. 6. 678820708-02 
CaTCh u" = atEs 2 .35313504E-03 
IMPULSe HFwI iPEMfcwTs 8,676928388 01 
TIME InTO FiJF.HTs 6,372V8775E 03 
DRHR3 6,44i35645 r J U 


coellihic u**‘t« i t placing vehicle in cor orbit 


cP3 ?.9»o33454"-02 
E P4 5 . 31=256993-04 
GAHA40 4 , < 7j 03849r-.il 


RPP4 6.85662935= 06 
T u aO 1 .62353487E 02 
VT4 7.61867509E 03 


R A P 4 6.86392630= 06 
R4 4.56375438E 06 
G A M T 4 O 9.24340470E-03 


AP4 6.86027/82E 06 
V4 7 .53710346E 03 
DELV4 9.70709611E 01 


the tpi irni’ion angle in relation to the taage t * 2 . 90000000E-01 
this stciroi' HfcTbRp.iNEs the catch up rate in the cuh orbit 

it ALSu SUnS UP the delta velocity AND catch up angle for the total mission 

TAUP4 5. 65 862208 ;3 PUPA40 6 . 366202 7.7E.02 DPDCUO 2 . 56925742F-04 DL 5R40 3.92219917E 00 
TTP I 2.367 06741:' .,4 O M R T O 3 . 66113665= ot D V J T 2 . 1 1 572 1 40 E 02 
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Appendix D 

PROGRAM MODULES AND DETAILED FLOWCHART 


D.l PROGRAM MODULES 


Name 

RK713 

RKG 

CONIC 

GMAT 


MAROT 

ARTAN 

POLY 

ECCV 


DEG 

FATT 

FATMU 

PRINT 


FATMUL 

TIME 

RANGA 

TRUE 


Function or Subroutine 

This is a seventh order Runge-Kutta integration routine which 
can integrate backward or forward 

The main routine of the integration package. Integration is 
variable step-size with an accuracy tolerance of 0.0000005 for 
the state 

Computes orbital parameters given the state. (Only for 
elliptical orbits.) 

Matrix transformation from space fixed inertial launch 
coordinate system (Xg) to the in-plane X' ' ' ' system 

X"" = C-iI-l t _0 N ^2 tA Z “ 2^1 *S 

Sets up elements of transformation matrix for an angle of 
rotation about the X, Y, and Z axis 

Arctangent from 0 to 2 tt or -ir to tt according to flag 
til 

Evaluates an n order polynomial given its coefficients 
Computes eccentricity vector e 



Earth's gravitational potential function. Evaluates the 
acceleration due to gravity for all three components 

Matrix transpose (3x3) 

Matrix multiplication (3x3 times 1x3) 

Calculates the U.T. in hours, min., sec, adds the U.T. to 
the mission Time "T", and prints out state and orbital 
parameters of each vehicle in flight. 

(Note: the program integrates in mission time, thus U.T. of 

launch is added to mission time from lift-off to obtain 
instantaneous U.T. time in flight) 

Matrix multiplication (3x3 times 3x3) 

Determines Keplerian time of flight between two positions on 
an elliptical orbit 

Computes range to and from descending node w.r.t. equator to 
the instantaneous radius vector 

Computes true anomaly from perigee to the instantaneous 
radius vector 
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4-1 1 TO PAGE 0-9 


RAP3 = R at - aHD + aHB, RPP3 = R cp , AP3 = 
(RPP3 + RAP3)/2. EP3 = ( RAP3 - RPP3)/ 
(RAP3 + RPP3) , RPP4 = R pT - aHD 

RAP4 = R at - aHD, AP4 = ( RPP4 + RAP4)/2 
EP4 = (RAP4 - RPP4)ARAP4 + RPP4) 


4 jjP3AP4(l - EP4 2 ) - EP4AP3(1 - EP3 2 ) 

R4 = AP4(1 - EP4 2 )/(1 + EP4 COS (e 4 )) 

V4 = j[2yAP3/R4) - y]/AP3j 1/2 

y 4 = TAN' 1 j(R4 EP3 SIN(e 4 )/(AP3(l - EP3 2 ))] 

VT4 = |[(2yAP4)/(R4) - y]/AP4| 1/2 ' 

V T4 = TAN' 1 |(R4 EP4 SIN 6 4 ]/[AP4(l - EP4 2 )]] 

AV4 = [VT4 2 + V4 2 - 2VT4 (V4) COS(y T/1 - yj] 1/2 
AVR = f (aHD) 14 4 


AP3 ( 1 - EP3 ) - AP4f 1 


FROM PG. D-7 


T.T 1 .X p .X p .X T ,X T Y 


T ] = TIME OF ORBITER INSERTION 
Rr P = |X p | . v vp = I I 
A p ~ ' V VP. R RP^ 

RrT = I I • ^VT = I ^T I 
a t = Rrt^C^ - V VT 2 R rp ) • 

ipAi -W/h f) 1 / 2 . * TA1 = [,/A t 3 ] 1/2 

A *1 = *PA] ' *tai 3 1/? 

A4> M ri = ^l^^P ; r 2“! 

H p = r rp , v vp cos . V E p = 1 ' aTIT 
y p = SIN' 1 t(X p • X p )/(R Rp • V vp )] P 

R AP ’ A P (1 + E P } * V AP = 


AP3 = (RPP3 + RAP3)/2, t p3 = 2n 

$PA3 = [A p 3 /y] 1/2 , a$ 3 = i pA3 - i TA1 , 
A* MR3 = A<^ o ( SFN02 ) t T 4 = T, + (SFN02) 
RAP3 = R AT - AHD + aSb. VPP3 = 

2~] 1/2 [2y RAP3/RPP3(RAP3 + RPP3] 1/2 , A V, = VPP3 


YES^ 

— ^Tlig^ ■ 






r h t 2 i 

■ E T ■ 11 ' A T y_ 


h t = R RT V vt cos y t 


r at = M 1 + e t^* t pi = 2w . ~ 

R PT = A y (l - E t ), AR T = R at - R pT 

T 2 = T 1 + T P1 /2 

JLIGHT = 0, SSFN03 = SFN03 
JULOS =_TULO 

x TS1 = x T , x psl = X p 
x TS i = x T , x psl = X p 


R CP 

= R AP’ v cp = 

Cp/R cp ] 1/2 

T CP 

= 2ir [R 3 p/y] 

1/2 

^CP 

= [y/R 3 p ] 1/2 



= *CP ' ^JAI ’ 

6 *MR2 

Ai 2 

(SFN01 ) T cp , 

AV 2 = 

V CP 

' V AP’ T 3 = T 

2 + 

(SFN01 ) T cp 
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4-3 FROM PG. D-10 


X T “ X TS1 


ORBIT INTEGRA- 
TION OF ORBITER 
AND SPACE STATION 
FROM T} TO 
TTEST1 



^AP ^RP l X pl ’ X * ^cpp ~ X p x X p 

V cp - X p x N cpp /( | X p x N cpp | ) , V cp = V Q V cp 


\ " v cp ’. X P’ v gl J V g I * x p " v cp ’ 

Xps = Xp, X ps = Xp, X TS = X T , X JS = X T 


, , / ORBIT INTEGRA- 

IM h 15 U TION OF ORBITER 

idj m SpACE STATI0N 

N/ \ TO TTEST2 

FROM PAGE D-1 9 ' — 1 
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LIFT-OFF TIME 
CORRECTION 


FROM PG. 0-11 


\ ^ 

ECCV 


/ — *■ 
-II 

e et 



FROM PG. D-1 3 


ORBIT INTEGRA- 
TION Or ORBITER 
AND SPACE 
STATION TO 
TTEST3 



TO PG. 0-13 


= |x t Iu/(2p - |x T r|x T | ). e t = |e et |. r atd = 

Aj(l ♦ Ej) - iHjj, Rpjg - A^(l - E^) - aH d 

E TD = !R AT0_' r pto ! { (r atd + “ptd* 

COSTTAPA = E et • (-X p )/[E T !x p | ] . R Ap * 

t 2R ATD R PTD^^ R AT0 4 R PTO^ 4 * R ATD ■ W C0STTAM ] 
♦ aHq, Rpp = | Xp | . £p n <R A p - Rpp)/(R A p + Rpp) 

Vpp * [ ( u/Rpp) ( 1 ♦ E p n ,/2 . N cpf - Xp x X p . 

Vpp = X p x N cp p/(|Xp x N C pp|). Vpp = Vpp Vpp 

V»pp- V V* |i g | - A P= (R PP^ R AP )/2 
.p, - 2n Ap 3/2 /t, TTEST3 = T ♦ t d1 • SFN02 



ORBIT INTEGRA- 
TION OF ORBITER 


6T = TG3 


AND SPACE 
STATION 


KPAS = 1 


PTA - PTASI 
TG3 


Tves 

PHTST x 0^ 


/\T0 PG. 

m o-u 


KPAS = 1 



aT = TG3 
JPAS ■ 1 


JPAS = 1 


ORBIT INTE- 
GRATION OF 
ORBITER AND 
SPACE STATION 

Ives 

-C TG3 > 10 > 
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D-14 
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LIFT-OFF TIME 
CORRECTION 
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4-17 FROM PG. D-18 


1 

r 

A4 e 
a* - 

A4 r 

1 

\ 



LAST PASS NODE CORRECTION 

aT L0 = A VV X TR1 

T X TS, 

T = TSTI, X T = X TS , 

X T = X TS 

TULO = TULO + aT L0 , 

aT LOT = 

TULOSS - TULO, TUTP 

= TUTP - 

aT L0 , TULOSS = TULO 



SPACE STATION ORBIT INTEGRATION (ONLY) 
FROM T = TSTI TO T = TSTI + aT, 

(IF aT LQ < 0 INTEGRATE 
BACKWARDS) 


'LO 


X TR2 " 

X T & T = TSTI + 

aT lo 

A *LOC = 

cos ' 1 [x TR1 • 

X TR2^I X TR1 1 l X TR2^ 

A4, loc = 

*- aT LO / 1 aT LO 1 

A4, loc 

aTTl = 

[A^oc ' (a *LOC 


' aT LO 

A* cu ]/ * 


TTEST2 

= TTEST2 + aTTl 




FROM PG. D-18 


TO PAGE D-18 


THIS IS THE SOLUTION 
RECYCLE AND PRINT 
OUT ISOLATED TIMELINE 
AND TARGETING PARAMETERS 



TO PAGE D-20 


ACCUMULATED VALUE* 

OF LAUNCH TIME CORRECTION 


X T = X TS’ 


*T " X TS* 


X p ” ^ps> ~ ^pg * 

T = TSTI 

*A t Lot = TULO - TULOS 


I SOLAS = 



I SOLAS + 1 




TO PG. D-l 1 


f MAJOR 
ISOLATII 
[ LOOP 
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Appendix E 

PROGRAM LISTING 


E-l 


. INC. 


M-1082 


° T I 


PROGRAM TARG 

PROGRAM TARG 

dimension xp(3>, xdp<3), xt(3>), xdt<3), xtri(3), xtr2(3)> tempk3> 

TBHP3(3 »' TE "»<<3>: T6MP5(3>, X0MEGX(3>. XAXU.s". BBB 

| <3 fL > i. cc E* 3 ' 3 ’; D ?®‘5: 3 ’; xpa <3 > » xts<3). xdts<3). sv<3> 

So XTSK3), XDTS1(3), XDP51 < 3 ) * XPSK3) 

DIMENSION A ( 7 ) 4 8(7), C (7). D<7>» E<7), F<7>, G<7>, H <7), Q(7), S( 

dimension xtlo(3), xdtlo(3) 

CFsl8S2 e 
Gn Qo 3W032E15 

GM2c 


PAGE 24 


s,72?gU585E»04 


ZEROsQ , 0 
ONE E l s 
7WOs2,0 

CwVa57, 29577951 
WA?OL8,1/Cmv 

KPASSsO 


* I SURF, I L I G 

- - -i A, 8, C, D, E, F, G, H, Q, S 

2100* t.tn.ts.ttol.hap.hpb* 

"»EN, XENON, THNN, ALFANiPHIN 

*es, xencs,thns» alfas»phjs 

s XLaMAL»BS\/D,A1*B1#Ci«T0LE,TY,DLHDiDL,MB»SFN0i,SFN02 


READ 2100 
DEAD 2100 
DEAD 2100 


UlNSsd 
DU8a§ 0 0 
DTs200 o 
Dt12s20oO 


JL l GHy a 0 
DySbDt 
PMIQbPHI 
AzOgAZ 
OELTHaDLHD 
OLHDsQI,HD*CF 

DUHSsQLHBdSF 

DMIsPhI/CNv 

AZsAZ/CNV 

XlAMAl>XLAMAl/CNV 

iSVDagSVD/CNV 

A1gA1/CNV 


XENGNa^ENCW/ONv 

7HNWsyMNN/CNv 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

U 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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PROGRAM TARG 

AlFANbALFAN/CNV 
PHINsPHIN/CNV 
XENCS«XENCS/CNV 
THNSsTHNS/CNV 
ALFAS»ALFAS/CNV 
PH I SrpH I S/CNV 
IF (TN-TS) 20 ,50*50 
20 IF <<fN-T>-TTOL) 40*30,30 
30 NORTHal 

tulo=tn 

. PRINT 1770 
GO TO 60 
A0 NORTHsO 
T U 1-0 a TS 

print i7eo 
GO TO 60 

50 IF < ( TS-T)-TTOL) 30, 40, 40 
60 IF (NORTH-1) 80,70,80 
70 TbTN 

at«an 

Et«EN 

XENCTbXENCN 
THNTafHNN 
ALFATbALFAN 
PhIIsPHIN 
GO TO 90 
80 TsTS 
At = AS 
Ef*Es 

xenct-xencs 

THNT=tHNs 
ALFATsALFAS 
PHI 1 = PH IS 

90 IF (JPASS-1) 100,140,140 
100 IF (THNT-PI/2,) 110,110,120 
110 AzL«Pi-ARSlN(COS(XENCT)/COS(PHi) ) 

AZO«AZL*CNy 
AZaAZL 
GO TO 130 

120 AZL«ARSIN<C0S (XENCT) /COS< PHI >) 

A Z*AZL 

azo«azl*cnv 

130 PRINT 1800 

phi io-phi i *cnv 
AlFATO=AlFA?*CnV 
XEnCtq«XEnCt*Cnv 
THNtO«tHNT*CNv 

PRINT 1790, T,AT,ET,xENCTO,f HN tO, ALFATO,PHI IO 
PRINT 1810, AzO 
140 PHITaPHI I-ALFAT 


PAGE 25 


A 

51 

A 

52 

A 

53 

A 

54 

A 

53 

A 

56 

A 

57 

A 

58 

A 

§9 

A 

60 

A 

61 

A 

62 

A 

63 

A 

64 

A 

65 

A 

66 

, A 

67 

A 

68 

A 

69 

A 

70 

A 

71 

A 

72 

A 

73 

A 

74 

A 

75 

A 

76 

A 

77 

A 

78 

A 

79 

A 

80 

A 

81 

A 

82 

A 

83 

A 

84 

A 

85 

A 

86 

A 

87 

A 

88 

A 

©9 

A 

90 

A 

9 1 

A 

92 

A 

93 

A 

94 

A 

95 

A 

96 

A 

97 

A 

98 

A 

99 


A 100 


E-3 
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T I D Y 


150 


170 


PROGRAM TARG 


call marqt <aaa,az-pi/ 2 .,i,i) 
call marot <bbb,phi, 3 ,«.d 
call mamul cccc,bbb,aaa) 
call marot <aaa,thnt, 2 ,-i ) 
call marot . (bbb,xenct#i,-d 
call mamul (Odd.bbb.aaa) 

Call mamul <aaa»ddd,ccc) 

Call marot c 9 bb»phit, 2 ,.i) 
call mamul <ccc,bbb,aaa > 

CALL FATT (DDD.CCC) 

KqmEGA(1)s§i N ( PH j J 

XoHEGA( 2 )a«O 0 S(PM! )«slN(AZ) 

XOMEGA( 3 )ssC 0 SCpHI )*COS{ AZ ) 

iF (KPASgei) 150,160,150 

ROsAT# $ ONE„ET*ET ) /(ONE*ET*COS(PHII)> 

VOsSQRT(GM* ( ’fwo/RO-ONE/AT) } 

C 0 N??WE T4N( ET * S,N ' PHI ‘ 1 ’ °«‘f=t,COS(PH!! , , 

TEMPI j[ 1 >«Ro 

?EMPl(2)aZER0 

TEMP1^3)»ZER0 

temps a >svo*sinigammo) 

TEMP2(2>*ZER0 

T|MP 2 ( 3 )oVO*COS(GAMMO) 

Call fatmu <xt,ddd,tempu 
gall fatmu cxdt»ddd«temp 2 ) 

BTYSssiQO.O 
MPASSJ bO 

( 5 r , (JPASS'II 170,370,170 
GALL RAnQA £ XT , X Dt , XQ MEG A ,p M If ) 

gall vCross ctempi»xy,xDt) 

Gall vCRQSg ( TEMP2, TEMPI. ynycflA. 
RRTaVMAG(XT) »XOMEGA) 

m ‘BARSlNCVDOT(XT,XQMEGA)/RRT) 

0«pSI*CNV * 

PiCDsCOSCPH! ) 

Pi£2)aSlN(PHi )*SIN(AZ) 

5 £ 3 )b«SIN(PMI )*COS( AZ) 


JRCOS ( TEMPS ) /*M«B ( TEMP S > ) 

_ ,-VD0 T <X0MEQA,TEMP4> 

BUMnDUMl*DUMl 

KgNC s^RTa N C SORT? ONE -BUM) ,DUM1 ’ a i 

PHlTOopHjT#CNV 

THNTOsTHnT^Oinv 

PR l MT 170 0 p pi i o 

iMPASSl-l) 180.300,180 

I) 200,190,190 
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A 101 
A 102 
A 103 
A 104 
A 105 
A 106 
A 107 
A 108 
A 109 
A 110 
A 111 
A 112 
A 113 
A 114 
A 115 
A 116 
A 117 
A 118 
A 119 
A 120 
A 121 
A 122 
A 123 
A 124 
A 125 
A 126 
A 127 
A 128 
A 129 
A 130 
A 131 
A 132 
A 133 
A 134 
A 135 
A 136 
A 137 
A 138 
A 139 
A 140 
A 141 
A 142 
A 143 
A 144 
A 145 
A 146 
A 147 
A 146 
A 149 
A 150 


Er -4 
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TIDY* 


PROGRAM TARG 


PAGE 27 


i90 NORTHrl 
GO TO 210 
200 NORTH.O 

southernly launch coefficients 

210 IF (NORTH-1) 220,230.230 
220 PSIDO«POLY(B # XENCO,6) 

PSIDsPSIDO/CNV 

DUM«C0S<XENC)/C0S(PSID) 

AA»SQRT(ONE-DUM*DUM) 

PHILS«aRTAN< <SIN(PSId)/COS(PS!d))#aA#-1> 

PhILS«PI2-PhILS 
dphRR*ph I LS-Ph i t 
GO ' T o 240 

NORTHERNS LAUNCH COEFFICIENTS 
230 PsIDo*POLY(A,XEnCo,6) 

PsIDspsIDO/CNv 

DuMaCOS<xENC)/COS(PSlD) 

AAaSQRTf ONE-DUM*DUM) 

PHILSbARTAN( ( SIN( PS ID)/C0S( PsId) > , AA# -1J 

PHILS.PI+PHILS 

DPHRRaPHILS-PHIT 

240 continue 

dphrro»dphrr*cnv 

PRINT 1710, PSIDO 
IF (PHIT-PI) 290.290,250 
250 IF (DPHRR) 29o,29o»260 
260 PHlDaSQRT(GM/(AT*AT*ATj ) 

TGR«DPHRR/PHID 
If (TQR-100.0) 270,270,290 
270 DTT3«10.0 

IF (TGR-20.0) 280.280,290 
280 DTT3af GR 
MP ASS I *1 
290 TFbt*DTT3 

caul rkg (phio, azo,xt#xDt,t,tf’) 

TbTF 

TuLODb 0 , 0 
PRINT 2210 

Call PRINT <TF,XT,XDT,AZ,PHj;TULOD) 

GO TO 170 
300 PRINT 1840 

PRINT 1820, PSIO 

PRINT 1830 , PSIDO 

XlDOoxENCO 

PRINT 1850, XENCO 

IF (ISURF-1) 310,330*310 

This is WHERE THE STEADY STATE SHOULD BE PROGRAMED IN THE. FUTURE 
310 DO 320 1*1,3' 

XP< J )*0 ,0 
320 XdPU ) *0 9 0 


A 131 
A 132 
A 153 
A 154 
A 155 
A 156 
A 157 
A 158 
A 159 
A 16 8 
A 161 
A 162 
A 163 
A 164 
A 165 
A 166 
A 167 
A 168 
A 169 
A 170 
A 171 
A 172 
A 173 
A 174 
A 175 
A 176 
A 177 
A 178 
A 179 
A 180 
A 181 
A 182 
A 183 
A 184 
A 185 
A 186 
A 187 
A 188 
A 189 
A 190 
A 191 
A 192 
A 193 
A 194 
A 195 
A 196 
A 197 
A 198 
A 199 
A 200 
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PAGE g 8 


PRINT I860 
Gq TO 400 

330 l F (NQRTH-1) 350,340,350 
340 AzNbPOLY ( C, x IDO, 6) 

Az Os &2 N 

THNNPsPOLY( D,X IDO, 6) 
TlN 3 PQLY(E,xiDO,6) 
PHIPNbPOLY( 0, X I DO, 6) 

PHI 1«PHIPN/CNV 

ALFATb 0,0 

THNT^THNNP/CNV 

AZ»AZN/CNV 

TiaTlN 

G0 TO 36Q 

350 AZSsPQLYIF , X I DO , 6 ) 

AZObA?S 

THNSPbP 0LY(G,XID0,6) 
TiSt3P0LY(H,XlD0,6) 
PHIPSbP0LY(S,XId0,6> 
phi i s ph i ps/cnv 

AlFATsO.O 

THNT8THNSP/CNV 

AZaA2S/CNV 

T % o T 3, § 

360 XENCT*X I OO/CN V 
PRINT 1870, A Z 0 

KPASIsi 


ROsNPgR*CF*RE 
RAqNAP&CF 
STb«RA'-ROI/{RA*RO) 
ATs-^RA^RO 5/TWO 
V0siQRT(GM$cON£*ET)/RO) 


370 IF CKPASS-lS 
3g0 DO 390 I b 1 , 3 
XPU )»XT(I ) 
390 XDPU S s Xq T ( i ) 
U 0 KP ASS a 0 


410,380,380 


PRINT I860, Xp , xDp 
PUNCH 2350, XP e XDp 
JpASSal 


@0 TO 60 
4?.0 TFsT^Tl 

00 420 1.1,3 
X710C I)«XT( I, 

420 XDTLOtl 5 »XDT<I ) 

PRINT 2320 

Cali rkg (phio,azo,xtiXDt,t , .tf> 


A 

201 

A 

202 

A 

203 

A 

204 

A 

205 

A 

206 

A 

207 

A 

208 

A 

209 

A 

210 

A 

211 

A 

212 

A 

213 

A 

214 

A 

219 

A 

216 

A 

217 

A 

218 

A 

219 

A 

220 

A 

221 

A 

222 

A 

223 

A 

224 

A 

225 

A 

226 

A 

227 

A 

228 

A 

229 

A 

230 

A 

231 

A 

232 

A 

233 


A 234 
A 235 
A 236 
A 237 
A 238 
A 239 
A 240 
A 241 
A 242 
A 243 
A 244 
A 245 
A 246 
A 247 
A 246 
A 249 
A 250 
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PROGRAM TARG 

CALL PRINT (TF,XT,XDT#AZ,PHr,TuL 0 D) 

PRINT 1840 
PRINT 1890, XT , XDT 

TULOSmTULO 
DO 430 1*1,3 

xtsi < i > *xt n ) 

XdTSIU >*XdT(I ) 

XPS1(! ) ^XP U ) 

430 XDPS1 ( | ) *XDP ( I ) 

440 RRP*SQRT(VD0T(XP,Xp >) 

WPaSQRT (VDoT (XDp.XDp ) ) 
ApsGM*RRp/(GM2oVVP*WP*RRP ) 
RRT*SQRT(VDOT<XT,XT)> 

WT»SQRT( VDOT(xDT,xDT) ) 

AT*GM*RRT/( GM2„VVT*VVT*RRT) 

PDPA1,SQRT( GM/< AP*AP* AP) ) 
PDTA1bSQRT(GM/(AT*AT*AT) ) 

DLPDIbPDPAI-PDTAI 

DLMR1«DLPD1*PI*SQRT<AP*AP*AP/GM> 

GAMAPbARSIN(VDOT(XP,XDP)/(RRP*WP)) 

HpnRRp*WP*COS(GAMAP) 

EPbSQRT(ONE-HP*HP/(GM*AP) ) 

RPPbAP* ( ONE-EP ) 

RaP»AP*(ONE*EP) 

VaP®SQRT(GM/AP)*80RT( (ONE-EP )/(ONE*EP>) 
GAMATbARSIN(VDOT(XT,XDT>/<RRT*VVT) > 
HT*RRT*VVT*COS(GAMAT) 

ET»SQRT(ONE-HT*hT/(QM*AT) ) 

RAT*At*<QNE*ET) 

RpT8AT*(ONE-ET> 

TAUPiap I2^SQRT< Ap*Ap*Ap/GM ) 

T2*Tl*T AuPl/TWO 

Drt*RAT-RPt 

GAMAPObGAMAP*CNv 

HaTb(RAT 0 RE)/CF 

PMDOTPbPDPA1*CNV 

HPTs(RPT»RE)/CF 

PHDOTTaPDTAl*CNV 

HAP*(RAP-RE)/CF 

DLPD10«DLPD1*CNV 

HPP*(RPP-RE)/CF 

GAMATObGAMAT*CNV 

DLMR10»DLMRi*CNV 

PRINT 2110 

PRINT 2120 $ RRP, VVP,GAMAPO,EP’,aP,HaP#HPP,PHDOTP 
PRINT 2130 

PRINT 2140, RRT,VVT,GAMATO»ET',aT ,HaT.MPT,PHDOTT 
PRINT 2150 

PRINT 2160, DLPD10,DLMR10iT2 


PAGE 29 


A 251 
A 292 
A 253 
A 254 
A 255 
A 256 
A 257 
A 258 
A 25? 
A 260 
A 261 
A 262 
A 263 
A 264 
A 265 
A 266 
A 267 
A 268 
A 269 
A 270 
A 271 
A 272 
A 273 
A 274 
A 275 
A 276 
A 277 
A 278 
A 279 
A 280 
A 281 
A 282 
A 283 
A 284 
A 285 
A 286 
A 287 
A 288 
A 2S9 
A 290 
A 291 
A 292 
A 293 
A 294 
A 295 
A 296 
A 297 
A 296 
A 299 
A 300 
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PROGRAM TARG 


RCP^RAP 

CRbSQRT<Gm/RCp) 

bSQRT<RCP*RCP*RCP/Gm) 

? AUCp op 1 2*DllM 

PDOTCpoONE/DuM 

DLPD2oPD0TCP-PDt Ai 

DLMR2oSFN01*DLPD2*TAuCP 

DBUV2aVCP«VAP 

T3sT2+SFN01*TAUCP 

PD0TC0bPD0TCP*CNVJDLMR20sDlMR2*CNV 
print 2170 

!F ! UL?o!i; ^: < V ";J^ CP ' PB 0teD,BLPD2,DL«R20,D E |.W.T3 

450 XNbZERO 
©0 TO 470 
460 KNnONE 
470 CONTINUE 

RPPSoRCP 

RaP33RPT-DLHD*DLHB 

AP3»(RPP3*RaP3)/TW0 

DUMoSQRT(AF3*aP3*aP3/GM) 

TaUP3bPI2*DUM 

PDPA3bONE/DUM 

DLPD3*PDPA3«PDTA1 

D^MR3bSFN02*DlpD3*TAuP3 

T4fflT3^8PN02*TAup3 

K3i;?iL s ^; RAp3/<ppp3 - <p ^ 3 * ppp ’"> 

print 2190 

PDPA30aPDPAS*eNv$DLPD30aDL.PD3*CNv 
0LMR30 S DLMR3 6 CNV 
30oDLM^3*CNV 

C!R e jLAR?’ATlorAT P C D DJ 3 A 0 L?t^ E 0 '- DELV3 - T4 ' DLMR3 ° 

RAP3BRAT-Df ! ,HQ*at.HB 
AP3 b£RPF 3*RAP3)/'TW0 
EP3s?RAP3®RPP3S/IRAP3'i'RPP3) 
nf?PT ®D|„HD 
®RAT-DtHD 
®(RPP4*RAP4 3/TWO 
iP4° ( RAP4 "RPP 4 ) / ( p a P4*RPP4 5 
P3PAP3*(ONg-gP3*gp 35 
P4 3 AP4*(ONg-gP4^gp 4 ) 

TH4aARS0SlCP3'»P43/(Ep3*p4 a cp4#p3jj 
R4aP4/(0f\jg*gp4*pQ5(TH4) > 

V4afORTI (QM2®Ap3/R4»QM)/AP3} 
@AMA4oARTAN(R4*gp3*SlN(TH4>'.pS’ .4 1 

VT4 q SQRTC(fiH2.A?4/R4iGN)/AM» ' l> 
eAMT4 B ARTA N IR4^Ep4*SlN(TH4), P 4;-l) 
fcUV 4 bSQRT(VT4*vT4*V4*V 4-TWO *VTA*V-4*C0S(GAMT 4-GAMA4 ) ) 


PAGE 30 


A 301 
A 302 
A 303 
A 304 
A 305 
A 306 
A 307 
A 30® 
A 309 
A 310 
A 3 U 
A 312 
A 313 
A 314 
A 315 
A 316 
A 317 
A 318 
A 319 
A 320 
A 321 
A 322 
A 323 
A 324 
A 325 
A 326 
A 327 
A 32Q 
A 329 
A 330 
A 331 
A 332 
A 333 
A 334 
A 335 
A 336 
A 337 
A 338 
A 339 
A 340 
A 34J 
A 342 
A 343 
A 344 
A 345 
A 346 
A 347 
A 348 - 
A 349 
A 350 
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PROGRAM TARG - 

TH40*tH4*CNV A 334 

GAMA40»GAMA4*CNV A 352 

GAMT40»GAMT4*CNV A 353 

PRINT 2210 a 354 

PRINT 2220, EP3,RPP4,RAP4, AP4,EP4,TH40,R4,V4 ( GAMA40 I VT4, QAMT40, DEL A 353 

1V4 A 356 

DuMbSGRT( AP4*AP4*AP4/GM) A 357 

TAUP4*P I 2*DUM A gg© 

PDPA4«0NE/DUM A 359 

PDPA40«PDPA4*CNV A 360 

DPDCU«PDRA4-PDTA'i a 361 

DPPHO«SLM*DELTH A 362 

DPHRsDPPHO/CNV A 363 

PRINT 2230 , DPPHO A 364 

DELVR*0 i 0 .a 365 

DLPMR4«DPDCU*TaUP4*<SFN03*XN)*dPHR A 366 

TTPlaT4*TAUP4*(SFN03+XNI*DPMR/DPDCU A 367 

DPDCUO«DPDCU*CNV A 36S 

DVI T»DEL V2*DELV3*DELV4*DELVR A 369 

DlMR4obDlPmR4*CnV A 370 

DpmRT1»DlMR1*DlmR2+DlMR3*DlPMR4 A 371 

DmRTO*«DPMRT1*CnV a 372 

PDP A4q»PDP A4*CNV A 373 

DPDCL)0«DPDCU*CNV a 374 

print 2240 * A 375 

print 2250 - A 376 

PRINT 2260, TAUP4,PDPA40,DPDCUO,DLMR40,TTPI,DMRTO,DVIT a 377 

call ranga (xp,xdp,xomega,phipi> a 37 s 

call ranga < xt, xdt, xomega, phitij a 379 

I 50LAS® 0 A 38§ 

I COR* 0 A 3§1 

PHlP10aPHlPl*CNV A 382 

PHIT10bPHIT1*CNV A 383 

PRINT 2270, PHIPIO.PHITIO A 384 

IF (PHIT1-PI) 480,490*490 A 365 

400 PHIT1«PHITi*PI2 A 3§6 

DpAl ®PH I Tl "■PH I Pi A 3@7 

GO TO 900 A 38S 

490 DPAlaPHlTl-PHlPl A 3©9 

500 IF (DPAl-OPMRTl) 510 , 520,526 A 3@e 

510 DPAl a DPAl*Pl2 A 3«i 

GO TO 530 A 392 

520 DPAlaDPAl A 393 

530 DTl«SFN03*TAUP4*TAUPl/TWO*TAUP3*SFN02 A 394 

DpHl»PLMRl*0LHR3*DpDCU*TAUP4*SFN03 A 395 

DpA10aDPAl*CNV A 396 

PRINT 1910, DpAlO A 397 

DpH2aDP A1-DPH1-DPHR A 398 

D T 2aDPH2/0LPD2 A 399 

TTESTlaTl*TAuPl/TwO-200, A 400 



^ 
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TTEST2aTl*TAUPi/TW0*DT2 

TSBST«T1*DT1*DT2 

print 1920 , sfnoi, sfno2, sfno3 

PRINT 2320 

CAU PRINT <T1,XT.XDT,AZ,PH!,TUL0> 

Caul RKG (PHlO,AZO,XT»XDT,Tl,fsSST) 

PRINT 2320 

CALL PRINT ( TSBST , XT, XDT, AZfPH! , TULO ) 

PRINT 2280 

DPHJObDPH1*CNV 

DPHt 0 sDPH 2 *CNV 

IF*(II ?r*J ! ^ Hi ? >DPH | 0 ' DT 2 - t S 8 ST,BT 1 ,Ti,TTSSTi,TTEST 2 

IF ULIG-1) 630*540,630 

TULQsiHUTTLE L I FT«OFF TIME IN SECONDS UNIVERSAL TtMf 
TV -number of days past JAN.i OF LAUNCH Year 
UMDAL ® LONGITUDE Of Launch Site 
INPUT CONSTANTS 

Omega b earths* rotation 

540 TbTSBsT 

550 DumsR I*XL A M A L”OMEGA*TULO 


70 


&QQ 


call marot (aaa.dum,2,®i> 

D(jMaAl*COS(Bl*Cl*TY) 

ALSV^DUM 

call marot (@bb,dum,3,=i> 
call mamul «ccc,bbb b aaaj 
call narot ?aaa p az»pi/two, i.d 
call marqt (bbb,phi,3,i) 
call mamul cddd.bbb, aaa> 
call mamul «aaa„ccc b dddi 
call fatt cb@b,aaas 

TEMPiiDoONE 
TEMPl(g)sZERQ 
TgMPl (3)aZER0 

CALL FATMU CTEMP2,BBB#TEMPi i 

UFLAGsO 

Call wgross (temps, xt,xdt) 

BSQoVMAQ ( T 6HP3 3/V0QT(XT,XT) 

call vuni? (?emp4»temp3) 

Call vcross c temps , tempa , temp2 » 
call vunj? (temps, temps) 

CALL VCROSS (TEMP2»TEMP5,TeMP4> 

BsVPTbARTAncVD0T(XT,TEMP5I,Vd0T <XT.TeMP2),1) 

IF CJLIGHT®!) 570 „ 1450 t 570 

VI f <@SypT«BsvD)- » 0001) 580,580,610 

IF (UFL a ® k 1) 590,628,390 

ChECKoBSVPT*BsD*DT 

If ( CHECK « 8svD ) 610,600,600 

UFLAQsl 

0?a(B§vD-@SvPT)/BSD 


A 401 
A 402 
A 403 
A 404 
A 405 
A 406 
A 407 
A 408 
A 409 
A 410 
A 411 
A 412 
A 413 
A 414 
A 415 
A 416 
A 417 
A 418 
A 419 
A 420 
A 421 
A 422 
A 423 
A 424 
A 425 
A 426 
A 427 
A 428 
A 429 
A 430 
A 431 
A 432 
A 433 
A 434 
A 435 
A 436 
A 437 
A 438 
A 439 
A 44 0 
A 441 
A 442 
A 443 
A 444 
A 445 
A 446 
A 447 
A 448 
A 449 
A 450 
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610 T2»T*DT 

CALL RKG <PHlO,AZO,XT,XDT,T, T2> 

T* T2 

GO TO 560 
620 PRINT 2320 

CALL PRINT <T2.XT.XPT,AZ,PHi’.TuL0> 

PRINT 2300 

PHSVOaPHSV*CNV 

ALSVOsALSV*CNV 

BSVPTObBSVPT*CNV 

AO*Ai *CNVSB03Bi*CNVSCObCi*CNV 

XLAMDObXLAMAL*CNV 

PRINT 2310/ AO.BO.CO.TY.XLAMDO'.PHSvO.ALSVO.BSVPTO 
TUTP«T 
DT3«T-TSbST 
DPHlaDPHl+DPDCU*DT3 
DPH2sdPa1-0PH1-dPMR 
DT2«DPH2/DLPQ2 
TTEST2bT1*TAUP1*DT2 
GO TO 640 
630 TUTPsTSBsT 
640 Do 650 1=1,3 
Xt<I)«XTSK!) 

650 XDT<I )«XDTSi< I > 

CALL rKG (PHIO, AzO,XT»XDT*Tl»fT E STl> 

CALL RKG (PHlO,AzO,xP*xDPrfTi.TTESTl> 

PRINT 2320 

call print <ttesti,xt.xdt, az.phi.tulo> 

PRINT 2330 

call print <ttesti,xp,xdp, az.phI.tulo> 

ToTTESTl 

660 CALL ECCV (GM.XP.XDP.TEMP3) 

Call TRUE (XP.XDP.TEMP3.PTA) 

Ep=VMAG(TEMP3) 

RpMaVMAG(XR) 

ApaRPM*GM/(GM2«VD0T(XDP.XDP)*RPM) 

CALL TIME (AP.EP.PTA.PI.GM.PI.TGPZ) 

IF ( TGP2 “ 3 0 • 0 ) 670 /670.680 
670 DT12*0NE 

IF (TGP2-TW0) 690*690.680 
680 T12 3 T+DT12 

call rkg <phio,azo,xt#xdt.t.ti2 > 

Call rkg <phio,azo,xp,xdp.t.ti2 > 

TaT12 
GO TQ 660 
690 RARbVmAGIXP) 
print 2320 

call print <ti 2 ,xt.xdt, az * phi V f ulo > 

Print 2330 

CALL PRINT (T12. XP . XDP » AZ .PHI »f ULO > 


A 451 
A 452 

A 453 

A 454 


A 497 

A 45© 
A 490 
A 460 
A 463, 
A 462 
A 463 
A 464 
A 469 
A 466 
A 467 
A 466 
A 469 
A 470 
A 471 
A 472 
A 473 
A 474 


A 476 
A 477 

A 47© 
A 479 


A 482 
A 483 

A 484 
A 485 
A 48© 
A 487 
A 48© 
A 4 
A 490 
A 491 
A 492 
A 493 
A 494 
A 499 
A 496 
A 497 
A 498 
A 499 
A 900 
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TST X»f 


PROGRAM TARG 


CAUL VGRQSS < TEMPI , XDP, XP > 
Call vcross <temp2,xp. tempi* 
call yunit { temp2, temp2) 

00 708 Ial,3 

780 TEMP3(I)bVC*TEMP2( I) 

00 718 j«l. # 3 

710 7IMR1C I)«TEMP3( !)-XDP( I ) 


PRINT 


PRINT 

CALL PRINT aSTI.XP.XDP.AZ.PHj.TULO) 

call ysRoss um.xp. temps* 
call eccv cgn s xp, temps, ev) 

PRINT 2018 c AM(i > ,AM( 2 ) # AMIS )»EV( 1 ) »EV( 2 ) »EV( 3 i .TEMPKD .TEMPI ( 2 ) , 


720 lot , 3 

i s' 

90 730 Ial fi 3 

i a ! 


XTStI )aXT<I 3 
730 KQYS { l ) sX0f C l 3 

740 CALL RK® <PHl0.AZ0#XT*XDT,T8T!’,fTEsT2) 
Call RKfi ^PHXO» AZO,XP.XDp,TSTi',TTEsT2) 

PRINT 2330 

prin/Imo <TTEs Ta>* p »x°P**z.PHi ,Tut-0> 
v a v«KI NT <'" f EST2,XT,xDT,4Z.PHi.T u l.O) 

To T i e§T c 

call eccv xt„ xdt, temp3) 

R?MqVMA@{ XT) 

A?sRTMo6M/C QN2-VD0T( XDT s XDT)*RTM) 

EtsVMAQ<TEMP3) 

RppsVMAG ( Xp 3 
RATDsaT* CONE*iT 3 °DLMD 
RP7QaAT*C0NE°ET )»DLMD 
ETDsCRATD^RPTDj/tRATD-s-RPTD) 

D0 750 I§1 0 3 
75 Q TgMP2 i 13 so^PCI 3 


CTYAPABVD0T«TgMP3.TEMP2)/CET*RPP) 

RaP s YH©^RAYO*RPTD/( <RATD*RPTD>*(RATD-RPTD)*CTTAPA >*DLHB 

gpa<RAP«RPP)/<R A P*RP P ) A * DLH B 

VPP s S@RT(C@M/RPP)MQNE*EP)) 
print 2330 

Call PRINT fTTEST2.XPiXDP.AZ.PHi.TuL0) 

PRINT 1930 mi. u t- 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


502 

503 

504 

505 

506 

507 

508 

509 

510 

511 
312 
513 

314 

315 

516 

517 

518 

519 

520 
321 

522 

523 

524 

525 

526 

527 

328 

329 

530 

531 

332 

333 

334 

335 

336 
537 

338 

339 

340 
541 

342 

343 

544 

545 

546 

547 

548 

549 

550 
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program targ 




CALL VCRQSS (TEMPI, XDP, XP ) 
CALL VCRQSS (TEmP3,XP,TEmP1) 
Call vunit (temp3,temp3> 

Do 760 1 = 1,3 

760 TEMP2( |)=VPP*TEMP3( ! ) 

Do 770 1*1,3 

770 TEMP1(I>*TEMP2(!)-xDP(I) 
VG2bvMAG(tEMP1) 

IF ( I SOL AS* 1 ) 790,780,780 

780 PRINT 2020 
PRINT 2000 
PRINT 2330 


call print <ttest2,xp, xdp, az.phi, tulo> 

CALL VCRQSS ( AM, XP, TEMP2) 

CALL ECCV (GM, XP, TEMP3,EV) 

PRINT 2010, AM(1),AM(2),AM(3),EV(1)#EV(2)»EV( 
iTEMPl ( 3 ) 

79o continue 

DO 800 1*1,3 
800 XDP ( l ) *TEMP2 ( I ) 

PRINT 1940 
PRINT 2330 

call print (ttests.xpixdp.az.phLtuLO) 
AP*(RPP*RAP)/TWO 


TaUP2»PI2*SQRT(aP*aP*AP/GM) 
TTEST3oT*(TaUP2*SFN02U25 0. 
JPASab 
KPASbO 


3) 


,TEMP1(1),TBMP1(2), 


Call RKG <PHIO.AZO.XT,XDT,TTESf2,TTEST3> 
CALL RKG (PHIO, AZO,XP,XDP»TTEST2,TTEST3> 
PRINT 2320 

call print <ttest3,xt,xdt, az.phi .tulq > 

PRINT 2330 

Call PRINT <TTEST3,XP»XDp,AZ*PHi.TuLQ> 

D T 22 = 5 , 0 
T =TT EST 3 

810 call ecc v (Gm,xp,xdp,tempd 
EPbvMAG(tEMPI) 

Call ECCV (GM,xT,xDT,TEMP4) 

Et»VMAG( TEMPA) 

RMTaVMAG( XP) 

Ap 8 RMT*GM/( GM2-VD0T( XDP, XDP)*RMT) 
RAP»AP*< ONE*EP) 

PpsAP* (ONE-EP*EP> 
DRTESfsRAP»VMAG(XT)*DLHD«>lOQ. 

CALL tRUE (XP, XDP, TEMPI, PTa> 

If (Et-TOLE) 820,940,940 
820 IF (DRTEST) 830,830,900 
830 IF ( P TA«P I ) 840,870,870 
840 CALL time ( AP»EP,PTA,Pf , GM,P!,TG3> 


A $ 5 % 

A 552 
A 593 
.A 994 
A 055 
A 596 
A 557 
A 99® 
A 999 
h 986 
A 961 
A 982 
A 963 
•.A 964 
A 969 
A 966 
A 967 
A 968 
A 969 
A 970 
A 971 
A 972 
A 973 
A 974 
A 979 
A 976 
A 977 
A 978 
A 979 
A 900 
A 981 
A 982 
A 9S3 
.A 984 
A 989 
A 986 
A 987 
A 988 
A 989 
A 998 
A 991 
A 992 
A 993 
A 994 
A 999 
A 996 
A 997 
A 998 
A 999 
A 600 


Et13 
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PRINT 1720, TG3 
IF (TG3-10,) 860,850*850 
850 Tl8sT+DT22 

call rkg <phio,azo,xp,xdp,t,ti8> 

Call RKG <PHI0,aZ 0,XT,XDT,T,T18> 
TsTl8 

IF UPAS-1) 810,870,810 

@60 DT22bTG 3 
JPASoi 
GO TO 850 
@70 RCPsVMAG CXP) 

PRINT 1730 
VCRcSqRT(Gm/RCp) 

call vcross (tempi, xdp, xp ) 

Call vcross <temp4,xpjtempi) 
call yuNir (tempa.tempa) 

DO 880 Iol,3 

880 TEMP2j I) bVCP*TEMP4< I) 

DO 890 Is I f 3 

890 TEMP3{ I)bTEMP2( I )-XDP( I ) 

VQ3n VMAG( TEMP3) 

GO TO 1140 

900 RS=VMAG(XT)-DLHD 

PTASIo-ARCOS( (PP”RS)/(EP*RS> )*PI2 
PRINT 1900 

call eccv <gm,xp,xdp,tempi> 
call true up, xdp, tempi, pta> 

PHTSTsPTASI-PTA 
IF (PHTST? @70,910,910 
910 Call time (AP,EP,PTA,PTASI, GM*, Pi , TG3 ) 
PRINT 1740, TG3 
IF <TG3-10. ) 930,920,920 
9gg Ti®§)T-j-DT22 

call rkg (pnio,azo,xt,xdt,t,t£8) 
call rkg <pmio,azo,xp,xdp,t,ti8) 

T@Ti® 

IF SKPAS-1) 810,870,810 
930 0T22at@3 
KPASal 

60 TO 920 
940 RmTbVmAG(XT) 

A?sR^t*Gm/(Gm2-VD0T (XD?,XDt )*RmT ) 

RATD«AT*<OnB*Et>-DlMD 

RPTD»AT*<0N6-Et)*DLHD 

A? Da ( rATD*RPT0> /TWO 

1t0b(RATD«3RPtD)/(rATD*RPtD) 

PYD®Af 0*< @NE»EtD*EtD) 

00 tso Ial,3 

950 ?^MP2(n° n TEMPl(I) 

eTAPAsVDQT^.TgMP2 5 TEMP4) / < ET*EP) 


A 601 
A 602 
A 603 
A 604 
A 605 
A 606 
A 607 
A 60S 
A 609 
A 610 
a 611 

A 612 

A 613 
A 614 
A 615 
A 616 
A 617 
A 618 
A 619 
A 620 
A 621 
A 622 
A 623 
A 624 
A 625 
A 626 
A 627 
A 621 
A 629 
A 630 
A 631 
A 632 
A 633 
A 634 
A 635 
A 636 
A 637 
A 638 
A 639 
A 64Q 
A 641 
A 642 
A 643 
A 644 
A 645 
A 646 
A 647 
A 648 
A 649 
A 650 
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RTDbPTD/( ONE*ETD#CTAPA) 
DRTESTaRAP-RTD- 100.0 
I PASS I aO 

IF (DRTEST) 1020,960,960 
960 Call VCROSS <TEMP2,XP,XDP) 

Call VCROSS <TEMP3,TEMP2,X0MEGA) 

call vcross (temps, temps, temps') 

S I NAPbVDOT(TEMPi, TEMPS )/<Ep*VM AG ( TEMPS ) ) 
COSAPeVDOT ( TEMPI, TEMP3)/(EP#VMAG( TEMPS) ) 
ALF APbARTAN(SINAP, COSAP,! ) 

call vcross <temp2,xt,xdt) 
call vcross <temp 3 ,temp2,xomeba) 

Call VCROSS <TEMP5,TEMP 3 ,TEMP2) 

si na To vdoT(temp4, temps )/<et*vmaG( temps > > 

COSAT8VDOT(TEMP4,TEMP3)/(ET*VMAG(TEMP3) ) 

ALFATbARTAN(SINaT,COSAT,1) 

DELAL«ALFAT-ALFAP 

DbETD*PP*C0S(DELAL )-PTD*EP 

Eb*ETD*PP*SIN(DElAL ) 

FspTD.pp 
Ab; ( E*E*D*D ) 

B = TW0 *F*E 
CcD*D-F *F 


RADs(B*B»4 i #a*C) 
sti = (.b,sqrt(Rad) ) / { two* a ) 
ST I Is ( -B+SQRTf RAD) )/< TWO* A) 
IF (STI) 970,970,980 
970 PTASIo-PI-ARSlN(STI) 


Go TO 990 

980 PTASIs-PI-ARSIN(STI I > 

990 Call TIME (AP,EP,PTA,PTASI ,GM’,pI,Tq 3) 
IF ( T G3sl0 , ) 1000,1000,1010 
1000 DT22=TG3 
IPASSIsl 
1010 Ti8sT*DT22 

call rkg (phio,azo,xp,xdp,t,ti8) 

Call RKG (PHIO, AZO, XT, XDT, T, Ti 8 > 

TsTj. 8 

IF < I PASS I -1 > 810 *1040,810 

1020 If (Pta*pd 1030,1030,1040 

1030 PTASIePI 
GO TO 990 
1040 RPMbVMAG(XP) 

PRINT 2320 

CALL PRINT (T18,XT,XDT, AZiPHl'.TULO) 

PRINT 2330 
Call PRINT 
Call ECCy (G 
RTM s VMAG(XT ) 

AtsRTM*GM/(GM2 


(T18,XP,XDP,AZ,PHI 
, XT, XDT, TEMPI) 


TULO > 


^VDOT(XDT,XDT)*RTM) 


A 651 
A 652 
A 653 


A 656 
A 657 


A 661 
A 662 
A 663 
A 66 4 
A 665 
A 666 
A 667 
A 668 
A 669 
A 670 
A 671 
A 672 
A 673 
A 674 
A 675 
A 676 
A 677 
A 678 
A 679 
A 680 
A 681 
A 682 
A 683 
A 684 
A 685 
A 686 
A 687 
A 688 
A 689 
A 690 
A 691 
A 692 
A 693 
A 694 
A 695 
A 696 
A 697 
A 696 
A 699 
A 700 
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EjaVMAG(TEMPl) 
rats A f*( one*Et ) 

RPTsAT*(ONE-ET) 

cau vcross < temp2,xdp,xp> 
call, vcross < temps, tempi, temp2> 

CTTAPPbVDOTc TEMPI, XP)/(RPM*ET) 
STTAPPbVD0T(TEMP3, XP)/(VMAG(TEMP3)*RPM) 

P?aAT* <ONE-ET*ET) 

RRTepf/(QNE*£T*CTTAPP) 

DELR1bRRT"RPM 
BeRPM«RAT“RPT 
IF (DELRl > 1050,1050,1060 
1050 BsoB 

1060 GsRAT*RPT*RPM/TWO*(CTTAPP*(Rpt = RAT)«*<RAT*RPT) ) 
Dlhisabs( (“B^sqrt (B*B”4,*C) ) /two ) 

DLH2»ABS( <-B-SQRT(B*B®4.*C) >/TWO) 
DELR«ABS(DELRi ) 

DRT1=ABS(DELR-DLH1) 

BRT2*aBS<DELR-DLH2) 

If (DRT1®DR72) 1070,1070,1080 

iqto dlhdpbdlhi 

GO TO 1090 
1080 DLHDPsDLH2 

1090 IF £ DELRl ) 1100,1100,1110 
1100 DLHDfa-DLHDP 
1110 RppsRPT-DLHDp 

RARsRAT°DLHDp 

Eps(RAR»RPp)/(RAp4.RPp) 

Aps£RAP*RPR)/TwO 

PPsAP*(OfyjE w EP*EP) 


1120 

115Q 

1140 

1150 


(GM/PP)* SORT (ONE* EP*EP* TWO* EP*CTTAPP) 
GAMMP e ARTA(\)( EP* ST TAPP, ONE+EP* CTT APP , - 1 ) 

Call vcross < tempi, xdp, xp> 
call vcross <temp2, xp, tempi) 
call vunit < temps, temp2j 
call vunit (tempi, xp> 

Do 1120 I«l,3 

TEMP2 ( n ^VR*COS ( GAMMP > *TEMP3 { I )*VP*SIN( GAMMP ) * TEMPI (I ) 

DO 1130 I a l a 3 

TEMPI < l ) sTEMP2? I )«XDP£ I ) 

V63aVMAG(TEMPl ) 

IF USQlAS«i) 1160,1150,1150 
PRINT 2030 


1 X NT 2330 

PRINT (Ti8,XP,XDP.AZ,PMl',TUL0) 

Call vcross um,xp,temp2 ) 
call eccv £gh,xp,temp2,ev> 

1TeMP 1£3) 10# AMfl)jAH(2),AM<3> ^V<l)#EV£2),EV(3),TEMPl<l),TEMPl(2), 


a 701 
A 702 
A 703 
A 704 
A 705 
A 706 
A 707 
A 708 
A 709 
A 710 
A 711 
A 712 
A 713 
A 714 
A 7i5 
A 716 
A 717 
A 718 
A 719 
A 720 
A 721 
A 722 
A 723 
A 724 
A 725 
A 7 26 
A 727 
A 72® 

A 729 
A 730 
A 731 
A 732 
A 733 
A 734 
A 735 
A 736 
A 737 
A 738 
A 739 
A 740 
A 741 
A 742 
A 743 
A 744 
A 745 
A 746 
A 7 47 
A 748 
A 7 49 
A 750 
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U6Q CONTINUE 

DO 1170 1=1,3 
1170 XDPU >®TEMp2(i ) 

I COR=I COR+1 
T«T18 

PRINT 2040 
PRINT 2320 

CALL PRINT <?18,XT»XDT,AZ,pHl‘,fuLO) 

print 2330 

Call PRINT <Tl8,xP,xDP,Az,PHr.TuL0) 

call rkg < phio, azo,xt,xDt, t, tutp> 
call rkg ( phio, azo, xp, xdp|t, tutP) 
print 2050 

PRINT 2320 

CALL PRINT (TUTP, XT.XDT, AZ, PHj'.fULO) 

PRINT 2330 

CALL PRINT (TUTP,XP,XDP,AZ,PHi',TUlO> 

IF (T18^TUTP*SSFN03*TAUP4-10 00 , ) 1200,1200,1180 
1180 SFN03sSFN03*,5 
DO 1190 1-1,3 
XT(I)«XTS1(I> 

XDT ( I ) aXDTSl ( I > 

XP( I )«XPS1 ( I ) 

1190 XDP<I)8XDPSl(I> 

PRINT 1950 
TULObTULOS 
DTaDTS 
GO TO 440 

1200 call vcross (tempi,xp,xdp> 
call vcross <temp2,xt,xdt> 
call vcross (temps, tempi, TEMP a’i 
watp»arcos (vdot( tempi *temp2)/(vmag (Tempi >*vmAq(temR2> ) > 
call ranqa (xp,xdp,xomega,phip> 

CALL RANGA (XT,XDT,XOMEGA,pHiT) 

call vcross (tempi. xt.tEmpsj 
call vcross <temp2, temps, tempi) 

SlNpNf ■VD0 t(TEMp2,XT)/(VMAG(tEMP2)*VMAG(XT) ) 

COsPNt«VDOt (TEMPS, XT )/(VMA-G(tEMP5)-*VMAG (XT-) ) 

PHINtsARtAN(SINPNt,COSPNt,1> 

call vcross (tempi, xp, tempsi 

call vcross (temps, temps, tempi j 

SINPNPbVDOTc TEMP2, XP)/( VMAG( TEMP2) * VMAG( XP) > 
CQSPNPsVDOT< TEMPS, XP ) / < VMAG( TEMPS }* VMAG( XP ) ) 
PHlNPaAR?AN<SlNPNP, COSPNP, 1) 

TEMP1(.1)*C08(PH! ) 

TEMPl(2)«SIN(PH!>*.SiN(AZ> 

TEMPI (3)b«S!N(PHI )*COS(AZ i 
call vcross <temp2, xt,xdt> 
call vcross (temps, temp2,xoMega) 

THNT»ARCOS(VDOT( TEMPI, TEMR3)/(VMAG( TEMPI ■)#VMAG( TEMR3 )■) ) 
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A 

m 

A 

752 

A 

753 

A 

7S4 

A 

793 

A 

756 

A 

757 

A 

79§ 

k 

759 

A 

768 

A 

761 

A 

762 

A 

763 

A 

764 

A 

769 

A 

766 

A 

767 

A 

768 

A 

769 

A 

770 

A 

771 

A 

772 

A 

773 

A 

774 

A 

775 

A 

776 

A 

777 

A 

778 

A 

779 

■ A 

78 0 

A 

781 

A 

782 

A 

783 

A 

784 

A 

785 

A 

7@6 

A 

787 

A 

788 

A 

789 

A 

7®0 

A 

791 

A 

792 

A 

793 

A 

794 

A 

795 

A 

796 

A 

797 

A 

798 

A 

799 

A 

800 
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1250 


1260 

1270 

12S0 


call vcross <temp2,xp,xdp) 

Call VCROSS CTEMP3,TEMP2,X0MEGa'> 

TnNPsARCOS< VDOT ( TEMPI » TEMP3 ) / ( VMAG ( TEMPI ) *VMAG ( TEMP3 ) ) ) 
DTHE»THNT*THNP ' 

DTHEOeDTHE*CNV 
PRINT 2060, DTHEO 
IF (WATP-WATOL) 1320,1210,1210 
1210 IF (PHINT-PI) 1220,1220,1270 
1220 IF (PHINP-PI ) 1240,1240,1230 
DELPHaPMlNT-PHINP 
GO TO 1420 

delph<»phint»phinp 

IF ( ABS<DElpH).pI ) 1260,1260,1250 
PHJNpapHlNp*p J2 
DELpHapHlNT«PHlNP 
GO TO 1420 
DELPH»PHIN T -PHINP 
GO TO 1420 

IF (PHINp-PIj 1290,1290,1280 
DELPHbPHINT-PHINP 
GO TO 1420 

delph.phint-phinp 
IF (ABS(DELPH) -PI) 1300,1300,1310 
DELPHbPHINT-PHINP 
GO TO 1420 

1310 PMINTbPHINT*PI2 

DELPHbPHINT^PHINP 
Go TO 1420 

1320 IF ( PH I T»P I ) 1330,1330,1370 
1330 OELPHbPHIT»PHIP 

If CPMIPePI) 1420,1420,1340 
1340 IF <ABS(DELPH)-PI) 1350,1360-1360 
1350 DELPHsPHIT-phip 
GO TO 1420 
PHITsPHIT+PI2 
DELPHsPHIT-PHIP 
GO TO 1420 

IF <PHIP»PJ ) 1380,1300,1410 
DELPHaPHI T-PHIP 

If i AqSIDELPH >«PI ) 1390,1400,1400 
0ELPHaPHI T -phip 
GO TO 1420 
PHI pspHlR*p j 2 

DELPHapHIT-PHIp 
GO TO 1420 

H10 DELPHapH IT -pH Ip 
1^20 DELPHObDELPH^CNv 
WATPObWAtPaCNV 
PHIP0aPHIP*CNv 
PHlTQsPHIT^CNV 


13?Q 

1300 

1390 


A 801 
A 802 
A 803 
A 804 
A 80S 
A 806 
A 807 
A 808 
A 809 
A 810 
A 811 
A 812 
A 813 
A 814 
A 815 
A 816 
A 817 
A 818 
A 819 

A 820 

A 821 
A 822 
A 823 
A 824 
A 825 
• A 826 
A 827 
A 828 
A 829 
A 830 
A 831 
A 832 
A 833 
A 834 
A 835 
A 836 
A 837 
A 838 
A 839 
A 840 
A 841 
A 842 
A 843 
A 844 
A 845 
A 846 
A 847 
A 848 
A 849 
A @50 
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1430 

1440 

1430 


1460 


1470 

1480 


1490 


1500 


1510 


PHINTObPH!NT*CNV 

PHINPObPH1NP*CNV 

PRINT 1970, WATPO, PHIPO, PH1TO, PH INTO, PHINPO 

PRINT 2070, DELPHO 

IF ULIG-1) 1460,1430,1460 

IF ( l SOUAS-1 ) 1460,1440,1460 

JLlGHTal 

GO TO 550 

BSVPAObBSVPT*CNV 

PRINT i960, BSVPAO 

GO TO 1650 

continue 

PRINT 2080, ICOR 

IF UCOR-2) 1470*1530*1560 

IF (ABS<DTHE)-.o2/CNV> 1530*1480*1480 

DLTLObDTHE/OMEGa 

TuIOsTUIO+DLTLO 

TaTST I 

PRINT 1980, DLTLO 
TuTPatUTP"Di.TLO 
DO 1490 I B 1 , 3 
XTRl ( i > “XTS < I ) 

XT < I ) aXTS < 1 > 

XDT < I >«XDTS < I ) 

DpHEEbDElPH-DpHr 

TDTL0«T*DLt1.0 

CAUL RKG (pHIO, AZO*XT*XDT*TSTi'.tDTLO) 

PRINT 2320 

call print <tdtlo,xT*xdt, az,phi',tulO) 

DO 1500 I»l,3 
XTR2( I)-XT< !> 


DLPLOC B ARCOS( VDOT( XTRl, XTR2) / f VMAGf XTR1)*VMAG< XTR2) ) ) 
DLPLOCb(DLTLO/ABS( DLTLO) )*DLPLOC 

DLTTl«(DLPLOC-( (DLPLOC*DPDCU) /DLPD2) -DLTLO* DPDCU)/DLPD2 
DLTT2«DPHEE/(DLPD2-DPDCU) 

DlTTT2aDlTTl«.DLTT2 

TTEST2aTTEST2*DLTTT2 

call ranga (xt,xdt,xomega,ph!tt) 
call vcrqss <tempi,xt»xdt) 
call vcross <temp2, tempi, xomega) 
call vunit ctemp4.tempd 

DuMiaVDOK XOMEGA, TEMP4) 

DuMaDUMl*DUMl 


X I NTbaRTaN ( 9QRT ( ONE-DUM ) , DUMt ,i ) 
TEMP1(1 > b COS(PHI > 

TeMP1(25sSIN(PhI)*SIN(aZ> 

TeMP1(3>b-SIN(PHI )*COS(AZ) 

ThTNT 8ARC os (VDOT« TEMPI, T eMP2)/(VMAG<TEMP2) > ) 

ThTNEbTHTNT^DThE 

CALL MAROT |AAA,A2-PI/2,,1,1) 


.A ©91 
A ©52 

A S3© 

A 854 


A '©3© 
A ©59 
■ A 

A ©61 
A 862 
A 863 

A 864 
A ©65 
A ©6© 
A 867 
A 868 
A 369 
A 870 
A 871 
A 072 
A 873 
A 874 
A 873 
A 876 


A ©79 
A §80 
A 831 
A ©32 
A 333 
A 884 
A 839 
A 33© 
A 887 
A 838 
A 8@9 
A 390 
A 891 
A 892 
A 893 
A 894 
A 095 
A 396 
A 897 
A 898 
A 899 
A 900 
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call marot (Bbb.phi , 3»n 
call mamul (Ccc.bbb, aaai 
call marot <aaa,thtne*2,-d 
call marot <bbb,x int,i.«i> 
call mamul <ddd.bbb, aaaj 
call mamul < aaa,ddd,ccc) 
call marot <bbb,phitt»2»*i> 

Call MAMgL (CCC.BBB, AAA) 

Call FATT (DDD,CCC) 

RRTaVMAG(XT) 

WTpVMAQ(XPT) 

GaMMA b ARS!N( VDOT( XT, XDT ) f < RRT* VVT ) ) 

TEMPI $ 1 ) sRRT 

TSMP1C2)bZER0 

TEMPI C3)sZER0 

call fatmu cxt.ddd. tempi) 

TEMPI (1)sVVT*SIN( GAMMA) 

tempi c 2 ) «zero 

TEMPl(3)aVVT*C0S(GAMMA) 

CALL FATMU (XDT. DDD, TEMPI) 

YpTST? 

DO 1520 1*1.3. 

XDP ( I ) sXDPS ( I ) 

XTS ( 1 5 = X T U ) 

XDTSC n=XDT( I ) 

1520 XPU )aXPSU ) 

TU“T«TULO 

IF (JINS) 1670.740. 1670 
1530 DPHEEbDELPH“DPHR 

IF (A8i^DPHEg)».05/CNV) 1580*. $580.1540 
1540 DTTT23DPHEE/CDLPD2-DPDCU) 

TsTST I 

Dq 1550 1*1,3 
XP ( I ) aXPS ( I 5 
XDP( I >«XDPS< I ) 

XT(I)*XTS(I) 

1550 XDT ( I 5 *sXDTS ( I ) 

TU®TST I*TULO 
TTESY2STTEST2+DTTT2 
GO TO 740 

1560 DPHEEbDELPH^DPHR 

IF { ABS(DPHEE)., 05/CNV) 1570.1570,1540 
1570 IF ( ABS(DTHE)- 0 02/CNV) 1580 , 1610 , 161 0 
1580 PRINT 2090 

IF USOLAS-l) 1590,1650.1650 
1590 ISOLASal 

DTULQfsTULO-TULOS 
PRINT 1980. DTULOT 
TeTST I 

DO 16QQ in fl 3 


A 901 
A 902 
A 903 
A 904 
A 90S 
A 906 
A 907 
A 908 
A 909 
A 910 
A 911 
A 912 
A 913 
A 914 
A 915 
A 916 
A 917 
A 918 
A 919 
A 920 
A 921 
A 922 
A 923 
A 924 
A 925 
A 926 
A 9 27 
A 92Q 
A 929 
A 930 
A 931 
A 932 
A 933 
A 934 
A 935 
A 936 
A 93 7 
A 938 
A 939 
A 940 
A 941 
A 942 
A 943 
A 944 
A 945 
A 946 
A 947 
A 948 
A 949 
A 950 
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P&66 4S 


XP(I)bXPSU) 

XDPU >«XDPS< I ) 

XT<I)aXTS(I> 

160 0 XDTU )bXDTS<I ) 

GO TO 740 

3.61Q DLTLObDTHE/OMEGA 
DO 1620 1*1,3 
XTR1U)«XTSU> 

XTU)bXTS(I) 

1620 XbT(Z >«XpTSCI > 

TULOafuiO*DLTLO 

TeTSTj 

PRINT 1980, DlTlO 
TUTP*tUTp-DLTLO 
Do 1630 I a l , 3 
XTR1U >*XTS< I ) 

XT C I ) «XTS 

1630 XDT < I ) «XDTS ( I ) 

Tr«TSTl*DlTLO 

call rkg (pmio, azo,xt»xDt«tst i'.tn 

PRINT 2320 

CALL PRINT (TF,xT,xDT, AZ.PHI.TULO) 

DO 1640 Ib1,3 
1640 XTR2( I)«XT( 1) 

DLPLOCbARCOS( VDOT( XTR1, XTR2)/< VMAG( XTR1)*VMAG( XTR2) ) ) 
DLPLQCb( DLTLO/ABS( DLTLO) )*DLPLOC 

DLTT1b(DLPLOC-( (DLPL0C*DPDCU)/DLPD2>-DLTL0*DPDCU)/DLPD2 
TTEST2bTTEST2*DLTT1 
GO TO 1510 
C 
C 

1630 TIMLAUbTULOS*DTULOT 
TIMINbTULOS*DTULOT*T1 

call rkg <phio, azo* xtlo, xdtlo, tulos# timlau) 

Do 1660 I B 1 # 3 

XT( I )bXTLO< I > 

1660 XDTU)bXDTLO(I) 

JiNS«t 

KINSbO 

DTHEbOMEGA*DTULOT 
GO TO 1510 

1670 If (KINS) 1680,1690,1680 
1680 PRINT 2340 

PUNCH 2350, XT.XOT 

punch 2350 , xt,xdt 

PUNCH 2100, XT, XDT 

Call PRINT m,XT»XDT,AZ,PHf»fuLoS+DTULOT> 

Go TO 1710 

1690 call rkg (phio,azo»xtsi,xdtsi»tulos*ti»tim!n) 

K X-NS 8 1 


A 931 
A 932 
A 953 
A 994 

A 935 
A 9S6 
A 997 



A 964 


A 936 
A 937 
A 968 
A 969 
A 970 
A 971 
A 972 
A 973 
A 974 
A 975 
A 976 
A 977 
A 97® 
A 979 
A 900 
A 981 
A 982 
A 903 
A 904 
A 985 
A 906 
A 907 
A 908 
A 989 
A 9.9 8 
A 99l 
A 992 
A 993 
A 994 
A 993 
A 996 

A 997 

A 998 
A 999 
AlOOO 
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1700 

ay 10 

1720 

1730 

1740 

1730 

1760 

1770 

1780 


Print 2370 
print 2360, xt.xDt 

CAUL PRINT <0.0, XT.XDT, AZ,PHl’,TtMLAU) 
DO 1700 1*1,3 
XT( I)«XTS1< I) 

XDT< I)«XDTSK I) 

Go TO 1510 
PRINT 2380, TIMLAU 


iaso 


i< 

ii 


FORMAT (1E16.8//) 

FORMAT (32H CIRCULARIZE W/O C D H EQUATIONS//) 

FQRMAt (1E16.8//) 

format (212) 
format <5Ei5.8/2Ei5.8> 

FORMAt (17H NORTHERLY LAUNCH//') 

format <i7h southerly launch//) 

format <?h timeeis.s.th atei5.8*7H etei5.b,.7h xenctoeis, 

18/7H ,THNTQE15.8,7H ALFAT0E15 , 6. 7H PHilOElS.B//) 

10 FORMAT (49H ORBITAL ELEMENTS OF SPaCE STATION IN-PLANE POINT/) 

.0 FORMAT <35H FIRST GUESS ON THE LAUNCH AZIMUTH*El5 .8// > 

1820 FORMAT <38H I NSTaNTANE I OUS LATITUDE OP INSERTION'S .8///) 

FORMAT <32H DESIRED LATITUDE POR I NSERT I 0N*E15 , 8/// > 

FORMAT «21H THIS IS THE SOLUTfoN) 

FORMAT <52H DESIRED VALUE OP INCLINATION POR TARGETING PURP0SE*E13 

1 , 8 //) 

60 FORMAT <61h INSERTION CONDjTiqnS DETERMINSD FROM STEADY STATE TRAj 

lECTORY//) 

7o Format M6h 

80 FqrhAt <40H 

1 0 7H YP£15.8,7M 
2 ZDPEi5,8//> 

90 format <59h state vector of space station 

1YION/7H XTE15.8,7H YTEl5,8,7H 
2H yDTE15 , 8 , 7H zDTEl5.8///i 
FORMAT < 39H INTERSECTION ASSUMMING CIRCULAR 0RB!TsEl5,8///> 

( 1®H PHASE ANGLE DPAlqElS. 8/ ) 

(3E16.8//) 

(33H BEFORE PERIGEE BURN AT TIME TTEST2/// ) 

(35H AFTER PERIGEE BURN AT T]ME TTEST2/// ) 

(72H1SFN03 HAS BEEN BUMPED BY ,5 BECAUSE THE TIME OF TP I OC 
BEFORE CDH///) 

(33H THE SOLAR VECTOR ANGLE ACH J E VEDbEi5 , 8// ) 

(7H WATP0E15.8.7H PHipOElS .8, 7H PH I TOE15 , 8 , 7H PH I NT0E15 , 

PH I NPE15 , 8// ) 

FORMAt <50H THE LAUNCH TIME SHOULD BE ADJUSTED BY DELTA TUL0*El5.8 

1 J 

1990 FORMAT (62H1TARGETING VALUES FOR THE CoV 100 NM C I RCULARZA T I ON AT 
1 APOGEE//) 

2 0 0 0 FORMAT C 29 H POSITION VECTOR FOR IGNITION//) 

§pio format I37H Targeting values for desired ellipse/24h angular momen 


ACTUAL LAUNCH AZIMUTH PROM ORBITEr !NSERTI0N»E15.8//) 

state variables of orbiter at Insert Ion/7h xpeis.s 

ZPE15.8/7H XDpEiS , 8, 7H vDpEl5,8,7H 


AT TIME OP 
ZTE15, 8/7H 


OrBITER INsEr 
XDTE19. 8, 7 


0 

1910 
1920 

1930 format 

1940 
1930 

ICURRED 

format 
1970 FORMAt 
l8,7H 


A 1 0 0 1 
A1002 
A1003 
A1004 
A1005 
A1006 
A1007 


A1009 

A1016 

Alois 

■Aioia 

Alois 

A1014 

Alois 

Alois 

A1017 

A1018 

A1019 

A1O20 

A 1021 
A1022 
A 1023 
A1024 
A102S 
A1026 


A102S 

A1029 
A 1 0 30 
A103S 
A 1 032 
A103S 
A1034 
A1035 
A1036 
A1037 
A1038 
A1039 
A1046 
A1041 
Al04f 
A1Q4S 
A1Q44 
Al04g 
A1046 
Al 0 4 7 
A1048 
A1Q49 
AlOiO 
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2020 

2030 

2040 

2030 

2060 

2070 

2080 

2090 

2100 

2110 

2120 


lTUM V(=CT 0 R/ 7 H AM (1 >El 5 .8 i7H aM< 2 )E 15 , 8 * 7 H AM( 3)613 . 8 / 20 H ECCENT 

2 RICJTY VECTOR/ 7 H E V ( 1 )E 15 , 8 . 7» E V < 2 >E 15 , 8 . 7 h EV( 3 )E 15 , 8 / 29 H VEL 

30 clTY TO BE GAINED VECT 0 R/ 7 H VG ( 1 )E 15 , 8 , 7 h VG < 2 >E 15 . 8 , 7 H VG( 3 )g 
419 , 6 //) 

FORMAT ( 42 H 1 TARGETJNG VAlUeS for The C 0 V PERIGEE BURN//) 

format < 46 hi targeting values for the cdh manuver For cov/> 

Format < 26 hicdh has been accomplished////) 

Format < 23 hiunive*sal time f 0 r TpI////) 

FoRMAf < 6 H DtHEoE 15 , 8 ///> 

( 7 H DELp H 0 E 15 , 8 // ) 

( 32 H 1 BEGIN NExT ISOLATION LOOP ICOR=l 12 ) 

( 21 H 1 THIS IS THE SOLUTION///) 

( 6 E 13 . 8 ) 

( 42 H PAPAMETERS FOR 30 x 10 O N,M, PHASING ORBIT/// ) 

< 7 H RRPE 15 . 8 . 7 H WPE 15 , 8 , 7 H GAMMAP 615 , 8 , 7 H EPE 15 , 

APE 15 . 8 . 7 H HAPE 15 , fl ' 4 7 H HPPElS. 8 , 7 H PHDOTPE 15 . 8 // ) 

( 32 H PARAMETERS FOR THE TARGET ORBIT///) 

< 7 H RRTE 15 . 8 , 7 H VVTE 13 . 8 , 7 H GAMAT 0 E 13 , 8 , 7 H ETE 15 , 

ATE 15 . 8 , 7 H HATE 15 . 8 , 7 H HPTE 15 , 8 , 7 H PHD 0 TTE 15 , 8 // ) 

( 74 H CATCH-UP RATE and ANGLE FOR THE HALP ORBIT OF! THE 50 X 1 


Format 
format 
format 
format 
format 
format 

18 / 7 H 

2130 format 

2140 FORMAt 
18 / 7 H 

2130 format 

100 NM PHASING ORBIT///) 

2160 format < 23 h orbital catch 

1 , 8 // 3 ?H TIME at APOGEE OF 

2 i 7 o format <sih first manuver 

2180 FORMAt < 7 H RCPEl 5 . 8 , 7 H 


UP RaTE»E 15 , 8 // 19 H ANGLE OF CATCH UPEE 15 
50 X 1 QO ORBIT»El 5 , 8 ) 

TO CIRCULARIZE 80 X 100 AT ITS APOGEE//) 

„ — , ..w. VCfElS 1 8» 7 H TAUCPElS ,8, 7 H PD 0 TC 0 E 15 , 

18/7H DUPD 2 E 15 , 8 # 7 H DLMRgOElS , 8 t 7H DELV 2 E 19 . 8 , 7 H TSE15.8///) 

2 i 9 o format (66h second burn transfer out of 100 nm circular towards 
1 COELLIPTIC///) 

2200 FORMAT < 16 H ORBITAL PER 1 0 D«E 15 , 8/I9H MEAN ORBITAL RaTE*E 13 . 8 / 15 H C 
1ATCH UP RATE-E 15 . 8 / 21 H IMPULSE REQU I REMBNT0EI5 » 8 / 18 H TIME! INTO FLI 
2 GHTPE 15 . 8 / 7 H DPHR 3 E 15 ,g/// ) 

2210 FORMAT M 5 m COeLLIPIc ORglT PL'aCINg VEHICLE IN CDH ORBIT///) 

2220 FORMAT ( 7 H EP 3 e 15 . 8 j 7 h RPP 4 e 15 . 8 * 7 H RAP 4 e 13 , 8 , 7 h AP 4 E 15 , 

18 / 7 H EP 4 E 15 . 8 , 7 H Th 40 e 15 , g, 7 H R 4 E 15 , 8 , 7 H V 4 e 15 . 8 / 7 H q 

ao ,« 2 i MA4 ?- 1 ^ 5 ,7H VT 4 E 15 , 8 , 7 H Q AMT 4 0 E 13 , 8, 7 h DEl V 4 E 15 . 8 //// ) 

2230 ^FoRMAT ( 50 h THE TP 1 I Qn I T I ON AnGlE IN RELATION TO ThE TAAGET*E 13,8 

2240 Format < 61 H this SECTION DETERMINES THE catch Up rate in the cdh 
1 - ORBIT/) 

2230 FORMAt ( 77 H It ALSO SUMS UP THE DELTA VELOCITY AND CATCH UP ANGLE 
1 FOR t HE TOTAL MISSION/) 

2260 FORMAT <7H TAgP4El3 t 8,7H P8PA40E15 . 8 , 7H DPDCgOElS, 8, 7H OLMR40E15 
18/7H TTPJE15, 8, 7H DMRT0E15, 8, 7H DVITEIB, §/// > 

2270 Format < 24 H RANGE ANGLE OF PURSUITaEl 3 , 8 // 23 H RANGE angle of TARGE 
1?bE15,8///) 

2280 FORMAt < 29 H COMPUTATIONS FOR SECTION 4 - 8 //) 

(7H DPH 10 E 15 , 8 5 7 H 0 PH 2 OE 15 . 8 , 7 H DT 2 E 15 , 8 , 7 H TSBSTE 15 , 

DT 1 E 15 . 8 , 7 H T 1 E 15 , 8 , 7 H TTEST 1 E 13 , 8, 7 H TTEST 2 E 15 , 8 /// ) 

( 5 QH COMPUTATIONS FOR LIGHTING CONDITIONS SECTION 4 - 10 //) 

< 7 H A 0 E 15 . 8 S 7 H B 0 E 1 S. 8 # 7 H C 0 E 13 , 8 , 7 H TYE 15 , 


2290 format 

18/7H 

2300 format 
2310 


A1051 
A 1052 
A 1093 

A 1854 

A 1095 

A 10 S 6 

A 10 P 7 


A 1662 

A 1 Q 63 

AH 


A 1066 

A 1067 

A 1068 

A 1069 

A 1070 
A1071 
A 1072 
A 1073 
A 1074 
Al0 79 
A1076 
A1077 
A 1078 
A 1079 
All 
All 
Aiosa 
A 1083 
A 1084 
A 1085 
A 1006 
A 1087 


A: 

A 1090 


A 1092 
A 1093 
A 1094 
A 1 Q 95 
A 1096 
A 1097 
A 1098 
A 1 0 99 
A 110 Q 
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2320 

2330 


10/7H LAMDA0E15.8,7H 


2350 


PHSV0E15 , 8» 7H AISV0E15 , 8 , 7H BSVPT0S15 , 8/ ) 
OF SPACE STATION//) 

OF ORBITER//) 

VECTOR OF TARGET AT INSERTION//) 


format <23h state 
format <i7h state 
Format <36h state 

FORMAt ( 6^13 1 6 ) 

FORMAT < 7 H XTE15.8,7H YTE15,8,7H 

l8,7H YDTE15.8.7H ZDTE15.8.//) 

70 FORMAT (35H STATE VECTOR OF TARGET AT LIFT OFF//) 

so format < 28 h the updated time of launchpeis . a//) 

END 


ZTE15,8,7H XDTE15. 


A110X 

A1102 

A1103 

AU04 
*1105 
A1106 
AU07 
All 06 
A 1109 

A1U0C 
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subroutine rk713 <to,tf,tol,xi.x,n,kt,m,8eta,alph,ch,ti,ab> 


SUBROUTINE RK713 (TQ,TF,TOL.Xi’.X,N,KT,M,BETA,ALPH,CH,TI, AB) 
SEVENTH ORDER RUNGE-KUTTA INTEGRATION WITH STEPSIZE CONTROL 
TF CAN BE GREATER THAN Tl OR LESS THAN TI AND RK7l3 WILL WORK 
M IS THE NUMBER OF STEPS NEEDED 

n is the number of differential equations 

Kt IS MAX NUMBER OF ITERATIONS 

ARRAY F STORES THE 13 EVALUATIONS OF THE DIFFERENTIAL EQUATIONS 
SUBSCRIPTS FOR ALPHA, BETA, AND CH ARE +1 GREATER THAN FEWLBERGS 
F(0) IN FEHLBERGS REPORT IS IN F < 1 , J > 

F ( I > IS IN F ( 1*1, J) 

FEHLBERGS REPORT REFERENCED IS NASA TR R-287 

parameters for deg subroutine must be stored in common 

DIMENSIONS MUST AGREE WITH NUMBER OF DIFFERENTIAL EQUATIONS AND 
NUMBER OF CONSTANTS IN THE PARTICULAR FEHLBERG FORMULA USED 
DIMENSION F ( 13 , 6 ) . XDUM ( 6 ) » TE(6), XI(6>, ALPH<13), BETA(13.12). 
1(6), CH ( 1 3 ) , AB ( 3 ) , ACCO ( 3 ) 

T = T 0 

DT=TF-T0 
M = 0 

DO 10 1=1, N 
10 X ( I ) a X I C I ) 

20 CALL DEQ (X.T.TE, AB,TI ) 

Do 30 1=1, N 
30 F ( 1 1 I ) aTE ( I ) 

DO 70 K = 2 * 1 3 
DO 40 1=1, N 
40 XDUM ( I ) =X ( I ) 

NNsK-1 
DO 50 I =1, N 
DO 50 J=1 , NN 

5 0 XDUM ( I )=XDUM( I ) + DT* BE T A ( K , J ) *F ( J , I ) 

TDUM=T*ALPH(K)*DT 

call deq <xdum,tdum,te,ab,ti> 

DO 60 1=1, N 
60 F ( K , I ) sTE ( I ) 

70 CONTINUE 
DO 80 I = 1 , N 
80 XDUM(I)=X(I> 

DO 90 1=1, N 
DO 90 L=l*13 

90 X( I ) = X( I )+DT*CH(L)*F(L, I) 

EPS=1. 

DO 120 1=1, N 

If all the variables being integrated have magnitudes whose 
absolute values are always much less THAN 1., then a value 
OF EPS LESS THAN ONE MAY NEED TO BE USED TO ACHIEVE AN ACCURACY 
AS SPECIFIED BY TOL. 

IF (ABS(XDUM(D)-FPS) 100,110,110 
100 AsEPS 

Go TO 120 


U 1 
C 2 
e 3 
C 4 
C 5 
Q 6 
G 7 
G 8 
C 9 
C 10 
G 11 
G 12 
G 13 
G 14 
X G 15 
G 16 
G 17 
C 18 
C 19 
G 20 
C 21 
C 22 
C 23 
G 24 
G 25 
g 26 
g 27 
g 28 
g 29 
g 30 
g 3i 
g 32 
C 33 
C 34 
g 35 
g 36 
g 37 
g 38 
g 39 
g 40 
g 4i 

£ 42 

C 43 
c 44 
C 45 
C 46 
g 47 
g 48 
g 49 
c 50 
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CTO,TF»TOL«XI,X,N#KT# 


m,beta,alph,ch,ti ,AB) 


110 A=XDUM(I) 

120 TEC I ) = D T ★ ( F ( 1 , I > + F ( 1 1 , I ) - F ( 1 ? , 
ER=ABS(TE(1J ) 1 

DO 140 152,1x1 

IF (ABS(Tf( I ) )-ER) 140,140,130 
130 Er= ABS ( TE ( I > ) 

140 continue 
dti«=dt 


I ) -F( 13 , I ) ) * 41 , / 84 Q . 


H = M*1 
AK=,8 

DTsAK*DT1*(T0L/ER)**,125 
IF (ER-TOL) 150,150,180 
150 T=T*DT1 

IF (ABS(DT)-ABS( TF-T) ) 170,17o’.160 
160 DTs TF - T A 0,1/0,160 

170 continue 
GO TC 200 
1 8 C DC 190 1=1, N 
loQ X( I )*XDUM< I ) 

20U IF (N-KT) 210,220,220 
210 IF (T-TF) 20,230,20 
220 TF=T 
230 RETURN 

End 


/A 


PAGE 8 


C 51 
C 52 
C 53 
C 54 
C 55 
C 56 
C 57 
C 58 
C 59 

C 60 

C 61 
C 62 
C 63 

C 64 
C 65 
C 66 
C 67 

C 68 

C 69 
C 70 
C 71 
C 72 
C 73 
C 74- 
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SUBROUTINE RKG ' rPHlL»AZ»XI»DXl,Tl,TF) ■■ 

SUBROUTINE RKG ( PH I L , A Z , X I , DX I *, T I , TF > D 1 

DIMENSION X<6), DX<6>, ALPH(1S)j BETA(13,12># CM(l3>. AB(3), XI(3) D 2 

j , DX I (3) D 3 

DO 10 1=1 #3 U 4 

X ( I ) = X I ( I ) b 5 

10 X( I +3 ) =DX I ( I ) 0 6 

GM=3 , 9660319E14 B 7 

RCONV=,1745329252E-01 D 8 

RPHIL*PHIL*RCONV 0 9 

HAZ=A7*RCONV 0 10 

AB(l)=SlN<RPHlL) 0 11 

AC3«C0S<RPHIL) 0 12 

AB(2) s-AC3*SlN(RAZ ) 0 13 

AB(3)»AC3*COS(RAZ) D 1=4 

DO 30 .1 =1 * 13 0 15 

DO 20 J=l,12 O 16 

20 BETA(I,J)b0, B 17 

■ALPH(I)=0, 0 18 

30 CH(I)sO, D 19 

CH(6)b34./105. 0 20 

CH(7)b9./35. O 21 

CH(8)«UH<7) O 22 

CH ( 9 ) =9 , /280 . D 23 

CH ( 1 0 } =CH < 9 ) . U 24 

CH(12)=4l,/840, 6 25 

CH ( 13 ) =CH ( 12 ) D 26 

ALPH<2)=2./27. B 27 

ALPH(3)=l./9, b 28 

ALPH(4)=l./6. O 29 

ALPH<5)=5,/12. b 30 

ALPH ( 6 ) = , 5 b 31 

AlPH(7)=5,/6. b 32 

ALPH(8)=l,/6. U 33 

A l. P H ( 9 ) = 2 , / 3 , D 34 

AlPH(10)=1./3. U 35 

AlPH(11>=1. D 36 

ALPH(13)b1, 6 37 

BETA(2,l)=2./27. D 38 

BET A ( 3 , 1 ) =1 . /36 .. D 39 

BETA(4,l)=l,/24. b 40 

BETA ( 5 # 1 ) =5 , /12 . D 41 

BETA(6,1)=.05 D 42 

BETA(7,l)=-25./108. O 43 

Bet a ( 8 , i ) = 3i . /3o o . b 44 

beta ( 9 , 1 ) = 2 , b 45 

BETM10, 1)b-91. /108, 0 46 

BETA(ll,l)s?383./4100. , D 47 

BETA(l2#l)s3./205. b 40 

BETA(13,1)=-1777./4100. 0 49 

BETA(3,2)=1./12. O 50 
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SUBROUTINE RKG (PHIl.AZ,XI»DXI,TI,TF) 



BFTA(4,3)=l./8. 

D 

51 

BETA(5,3)=-25./16. 

D 

52 

BETA(5, 4 )=“BETA(5,3) 

U 

53 

beta <6* 4 >s,;»5 

U 

54 

BF:TA(7>4)=125./108. 

U 

55 

BFTA<9,4)=-33./6. 

U 

56 

BFTA(10,4)s23./108. 

L) 

57 

BETACU, 4 )5-341. /164. 

U 

58 

BETa(13 # 4)=8ETa(11, 4 ) 

0 

59 

bet a ( 6 « ^ « 2 

0 

60 

BFTA(7 P 5)2-6b„ /27„ 

L 

61 

BFTA(8, 5>)a61. /225. 

D 

62 

BETA(9, 5)=704„/45. 

D 

63 

BETA(10,5)s™976./135, 

D 

64 

BETA<11,5)»4496./1025. 

U 

65 

BETA(l3,5)«BtTA(ll,5) 

u 

66 

BETA(7,6)ai2b./54. 

u 

67 

BETA(8#6>*-2./9. 

6 

68 

BETA(9, 6)3.107. /9. 

0 

69 

BFTA(lO»6)s311 t /54. 

u 

70 

6tTA(ll, 6)3-301. /82 , 

u 

71 

BFT A ( 12 , 6)3-6. /41. 

0 

72 

QETA(i3 } 6)s^2H9„/82, 

D 

73 

BFTA(8,7)s 13 o /900 . 

0 

7 4 

BETA(9 s 7)s67./90. 

D 

75 

3ETA(10,7)3-19 o /60 . 

0 

76 

3ETA(il P 7)s2133./41!30. 

D 

77 

BETA(l2p7)a-3./205. 

D 

70 

BFTA(i3p7)s2193./4100. 

t) 

79 

BETA(9, 8)s3, 

0 

@0 

3ETA(lO»8)sl7„/6„ 

D 

81 

BETA(H»8)a4S./82. 

y 

82 

BETA(12, 8)3-3. /41. 

c 

83 

8ETA(13»8)a51./82. 

u 

84 

BETAdO,9>8-l./12. 

0 

85 

3ETAU1, 9 ) 345. /164. 

D 

86 

BETA(i2,9)b3./41. 

U 

87 

BeTA(3.3c9)s33,/164. 

0 

80 

BETAdl, 10)318, /41. 

U 

89 

BETA(l2ilO)«6./41. 

u 

90 

BETA(13»10)b12./41. 

u 

91 

BtTA(l3il2)al. 

u 

92 

CALL DEQ (XfTI ,DX,A9,TI) 

D 

93 

ToL-,bE°06 

D 

94 

T 0 - T I 

D 

95 

CALL R K 7 1 3 (TO,TF,TOL#XiX»6i20nO,M,BETA»ALPH*CH,TI#AB) 

D 

96 

call nta tx. rF t nx» ab, ti j 

U 

97 

DO 40 1=1,3 

U 

90 

X I ( I )sx ( I > 

u 

99 

4 0 UxI(I)sX(I+3) 

u 

100 
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SUBROUTINE RKG -(PH1L»AZ,XI»DXI,TI,TF) 

RETURN 0 101 

End d 102 ° 


E-29 



NORTHROP SERVICES. INC. 


M-1082 


Subroutine conic <r, v,az.phi,aa,ap,enc,thtn,th,e,p,a,alfad.ra.rp,c PA - E 
,-< U ou?; J ! INE CUNIC • r ' v ' aZ ' PhI * aa *ap, e nc,thtn,th,e,p,a,alfad,ra,rp,c 

!o » PH I I ) 

1UU (^3 ) ^PUf 3) ^yi 3) p H<3> ' V ! 3> ' HU<3)# t HNV(3), THNU<3>. 

1UU(^), "U < 3 ) » XI (3), T ( 3 ) , SIJ(3), B(3), CU(3) 

W( 1 )=SIN(PHI, ) 

W ( 2 ) s-COS (PH I ) *S IN ( A 2 ) 

W ( 3 ) =CQS (PHI )*CoS(AZ) 

Call vunit (ru,r> 

call vcross <h,r, V ) 

call vunit <hu,h) 

Call vcross <thnv,h,w) 

call vunit <thnu,thnv> 

call vcross (ou.hu, ru) 

call vcross <pu,thnu,hu> 

GH=3.93603iP7PEi4 
RM = SQRT(VDOT(R,R) ) 

PsVDOT<H,H)/Qm 
Rp=VD0T ( V # RU ) 

AsGM*RH/(2,*GM-RM*VD0T(V,V)) 

TESTs( 1 .-P/A> 

IF (TEST) 20 , 20,10 
10 EsSQRT(TEST) 

GO TO 30 
20 Es 0 , 0 

30 Continue 

C 0 STHs(P«RM)/(E*rM) 

SINTHa(RD/E)*SQRT(P/GM) 

DO 40 1=1,3 

XI( I ) s RU ( I)*COSTH=QU( I ) * S I NTH 
ALFAD=ARTAN( VDUT( XI,PU>, VD0T(XI, thnu>,d 
THsARTAN(SINTH # COSTH, 1 ) 
call vcross <cu,hu, thnij> 

Ta)=cos(PHir T<CU ' RU) ' VDOT<THNU,Ru) ' 1> 

T(2)«SIN<AZ)*SIN(PHI ) 

T(3)=-C0S(AZ)*SIN(PHI) 
call vcross <su,w,n 

rl , | T | N 3 f?Dici VOOT(THNU ' Su> ' VDOT(THNU,T),D 

Call vcross <r,w,thnu) 

Re = 6378U6!° 0T<HU ’ S ’''' DOTlHU, “’'' 1 ’ 

Cwv= 1852 . 

C3s*Gm/A 
RAsA*( 1 ,+E) 

Rp=a*<i,.E) 

AAs(RA-RE)/CNV 
ARs(RP-RE)/CNV 
return 
end 
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40 


fc 

t 

h 

fc 

t 

t 

b 

t 

fc 

e 

fc 

fc 

fc 

i 

£ 

fc 

fc 

fc 

fc 

fc 

fc 

fc 

fc 

fc 

E 

fc 

fc 

b 

fc 

fc 

fc 

fc 

fc 

fc 

fc 

fc 

b 

fc 

fc 

£ 

fc 

fc 

fc 

£ 

fc 

fc 

b 

t 

fc 


3 

4 

5 

6 
7 
6 
9 

10 

u 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 ' 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


44 

45 

46 

47 

48 
49 - 
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SUBROUTINE GMAT ( PH I . AZ I , THTN . ENC , G , P I ) 


SUBROUTINE GMAT ( PH I , AZ I » THTN, ENC# G » P I > 

D<3 - 3 >' TE(3 - 3 >' SC3.3J 

call marot (B,phi,3,i> 
call marot cc, thtn, 2, »i > 
call marot <d,enc,i,-d 
call mamul <g,b,aa) 
call mamul <te,c,g> 

CALL MAMUL (G,D.TE) 

RETURN 
E ND 


M-1082 
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F 1 
F 2 
F 3 
F 4 
F 5 
F 6 
F 7 
F 8 
F 9 

F 10 
F 11 ® 
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"TIDY* 

SUBROUTINE MAROT ( A , ANGLE , K , L ) PA 9 E 22 

subroutine marot ( a, angle, k,l> 

DIMENSION A (3,3), C<3»3) * ^ 

SaNG=S I N (ANGLE ) H 2 

CANGsCOS(AngLF) H 3 

DO 10 1=1,3 H 4 

DO 10 J = 1 , 3 H 5 

10 A ( 1 f J ) sO , Hg 

M=(-3wK**2+11*k=4)/2 H 7 

N s ( 3*K**2-I3*K+3 6) /2 R 8 

A ( K, K ) a 1 , 0 H 9 

A(M 8 M )s CANG H 10 

A(N 5 N)sCANG H 11 

A(M,N)aSANG H 12 

A(N,M)a-SANR h 13 

IF <L> 20,20,50 H 14 

20 DO 30 1=1,3 N 15 

DO 30 J=l,3 H 16 

3ii C(I,j')aA<J,n H 17 

DO 4Q 1 = 1,3 H 18 

DO 40 J = 1,3 H 19 

4 II A ( I , j ) =C ( I , J ) H 20 

50 Return h 21 . 

END h 22 

H 23- 
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* T ID Y * 


FUNCTION ARTAN { SANG , GANG, I SW > 


C 

C 

c 

c 

c 

c 

c 


function artan (sang.cang, iswj 

this subroutine use the sine and cosine of the function 

and places the angle in the proper quadrant, 

IF ISWbI the angle is put BETWEEN 0 AND 2 PI 
IF JSWB-1 THE angle IS PUT BEf WEEN - PJ AND + Pi 

Pls'3, 14159265 
IF (SANG) 1,7,10 

1 IF (CANG) 2,3,4 

2 ArtANs-P I+ATAN(SAnG/CAnG) 

GO TO 5 

3 ArtANs-P 1/2, 

Go TO 5 

4 ARTANsATAN (SA nG/CAnG) 

5 IF ( I s W > 14,14,6 

6 ARTANs2,*Pi*ARTAN 
GO TO 14 

7 IF (CANG) 8,9,9 

8 ARTANsPI 
GO TO 14 

9 ARTANbQ, 

GO TO 14 

10 IF (CANG) 11,12,13 

11 ARTANsRI+ATANCSANG/CaNG) 

GO TO 14 

12 ARTANsPI/2, 

GO TO 14 

13 ARTANaATAN ( SANG/CANG ) 

14 RETURN 


end 


(Blank card) 
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L 1 
L 3 

L 4 
L 5 
L 6 
L 7 
L 8 
L 9 
L 10 
L 11 
L 12 
L 13 
L 14 
L 15 
L 16 
L 17 
L 18 
L 19 
L 20 
L 21 
L 22 
L 23 
L 24 
L 25 
L 26 
L 27 
l, 26 
L 29 
L 30 
L 31 

l 32o 
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M-1082 


* T I D Y * 

FUNCTION poly (CiX»N) 

function poly <c,x,n> 
c c is the coefficient array 

c x is the independent variable 

C N IS THE D E 15 R E E OF THE POLYNOMIAL 

DIMENSION CU> 

POLYSO .0 
KsN*l 

10 POLYsC (K )+P()LY*x 
K = K®1 

If c k 0 g t o > 10,20 

20 RETURN 
End 


PAGE 2 © 


■-J 1 

'j 2 

: j 3 

<j A 

! j 5 

j 6 

j 7 

j 3 

j 9 

'U 10 

j U 

1 j 12o 
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3200 FORTRAN (3,0)/RTS 


360 


SUBROUTINE ECCV(GM,XP,XDP, TEMPI) ‘ „ 

DIMENSION XP(3)»XDP(3),TEMP1(3) , TeMP 2 < 3 ) , TEMP3 < 3 > 
CALL VCROSS (TEMPI, xP, XDP) 

CALL VUNIT <TEMP?,XP> 

CALL VCROSS (TEMP3, TEMPI, XDP) 

DO 360 1=1,3 

TEMRK I )a-(TEMP2( I )*TEMP3( I >/GM> 

RETURN 

END 


407 

40® 

409 

410 

411 


FORTRAN DIAGNOSTIC RESULTS FOR ECCV 
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"" I 1 


tidy* 


SUBROUTINE DEQ (X,T,DX, AB,TI ) 


PAUE 


subroutine peo <x,t,dx, ab.ti ) 

DIMENSION X ( 6 ) , D X ( 6 ) , AB<3} # XDUM(6), ACC0(3) 

GM = 3,9B6031<?E14 
AA=,6378l66E*07 
Fj=l,62345E-03 
FH=-5 . 75E-06 
FD=7. 0756-06 
DO 20 1=1,3 
20 Dx( 1 )»X( 1+3) 

R2=X(1)*X(1)+X(2)*X(2)+X(3)*X(3) 

RsSORT (R2) 

RI=1./R 
R2I*1./R2 
B=AA*AA*R2l 
BBsA A<»R I 

A«(AB(1>*X(1)+AB(2)*X(2)*AB(3)*X(3))*RI 

A2s A*A 

A4sA2*A2 

GRaB*(FJ*U.-5.*A2)+3.*FD*(i./7.-2.*A2 + 3.*A4)*B4-FH*BB*A*(3.-7.*A2) 

1 > 

GP = B*(2,*Fj + a + 4,*FD*A*(3./7,oA?)*B + 3,*FH*BB*(A2-1./5. ) ) 

DO 30 1=1,3 

30 Dx< I+3)«-GM*R2I*< <1,+GR)*X( I )*R!+GP*AB< ! > > 

RETURN 

End 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


2 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
25° 
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subroutine fait (bbb*aaa) 


subroutine fatt <rbr,aaa> 

DIMENSION BHH ( 3 # 3 ) , AAA(3,3) 
DO 10 l. 3 1 * 3 
DO 10 M*l#3 
10 H B b ( L , M ) = 0 , 

DO 20 J = 1 # 3 
DO 20 1=1,3 
20 BB8 ( J , l ) ■ AAA l I , J) 

RETURN 

END 
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G 1 
G 2 
G 3 
G 4 
ci 5 
6 6 
t i 7 

G 8 
G 9 
G 10® 
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TIDY* PAliE 28 

SUBROUTINE FATMU (EEE, AAA.DDD) 

subroutine fatmu (EEE» aaa.ddd) r i 

DIMENSION fcEt ( 3 ) , AAA(3,3), DDD < 3 ) R 2 

DO 1C L = 1 » 3 R 3 

10 E E t ( L } b 0 , R A 

DO 20 1=1,3 R 5 

DO 20 J= 1,3 K 6 

20 EEE(I )bEEE(D+aaa( I, J)*DDD(J) R 7 

RETURN r 8 

End k 9 ° 
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* T I D Y * 


SUBROUTINE PRINT <T,RI,VI,AZ,PH,TULO) 


PAUE 


SUBROUTINE PRINT ( T , R I , V I , AZ , PH , TULO > 

DIMENSION RK3), V I ( 3 ) 

TULOIoTULO*T 
TT = T 
I C 0 R a o 

10 HRxTT/3600. 

IHRbHR 

XMINs(H-IHR*3600 . )/60 . 

HlNsXMIN 

SECaTT- I HR* 36 0 0 . -MIN*60 , 

IP (ICOR-1) 20,30,30 
20 PRINT 40, I HR , M I N , SEC 
TTsTULOI 
I COR = l 
GO TO 10 

3C PRINT 50, IHR,MIN,SEC 
CnV«57. 295779513 

call conic <ri , vi, az,ph,aa. ap',enc.thn»th,e,p, a,alf#ra,rp#C3,pwi i > 

fcNCl=ENC*CNV 

THN1=tHN*CNV 

TH1»TH*CNV 
ALFl*ALF*CNV 
PHI IlePHI UCNV 

print 60, t,ri , vi , aa, ap,ra,rp,p, a,e*C3,enci,thni,thi, alfi.phi u 

RETURN 


4o format (2X,ivh time from lift-off/sh hrs*, I2,3x,5h min=, i2,3x,5h s 

lEc*i El5 , 8 // > 

50 Format ( 2 X , 1 5 H UNIVERSAL T I ME/5H HRS = , I 2 , 3X , 5H M I N* , I 2 , 3 X , 5H SEC* , 
1E15.8//) 

60 Format </3x,4HTIMF.Ei5.8/5x, 2H X,El5,8,6X*lHY,El5.8,6X,HZ,El5.8,5 
1X,2HXD,E15.8,5X,2HYD,F15.8,5X,2HZD,E15.8/4X,3H A A , El5 . 8 , 4 X , 3H AP,E 
215.B,4X»3H RA.E15.8,4X,3H RP . El3 . 8 . 5X , 2H P , E15 , 8 * 2X , 5H A.E15.8, 
3/5X,2H E,Ei5.8,4X,3H C3,E15,3,4X,3 hENC,Ei5.8,4X,3HTHN,Ei5.8,5X,2HT 
4 M , E 1 5 . 8 » 2 X , 5 H A L. F A D , E 1 5 . 8 / 7 H PHII0E15.8///) 
end 
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TIDY* PAbE 30 

SUBROUTINE FaTML (CCCtBBB, AAA) 

SUBROUTINE FATML ( COC > RB8 » A A A ) L 1 

DIMENSION CCC(3,3), BBB(3.3>. A A A < 3 ' 3 } L 2 

no to L=l»3 I 3 

DO 10 h = i>3 l 4 

1 0 C C C ( L . M ) s 0 , L 5 

DO 2 C J = l,3 L ^ 

DO 20 1=1,3 l 7 

DO 20 K=l,3 L 8 

20 CCC( i , J ) s C C C ( I , J ) + B H B ( I,K)*AAA(K,J) L 9 

NET U H U L 10 

Emu l 11° 
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tidy* 


SUBROUTINE TIME ( A , E , THA , THB , OM , P I , TF ) 


SUBROUTINE TIME ( A , E , THA , THB . QM , P I , TF ) 


THIS SUBROUTINE DETERMINES THE KEPLERlAN TIME 

between two positions on an ELl'PTlCAL orbit 


OF FLIGHT 


10 


20 
3 0 


DIMENSION TH ( 2 ) * S I NE C 2 ) « 
TH(l)rTHA 


COSE ( 2 ) , EC A ( 2 ) » X M ( 2 ) 


TH ( 2 ) sTHB 

00 10 1 = 1,2 


SlNE(I)sSORT(l,-E**2)*SIN(TH(i))/(l.*E*COS<'TH(I))) 
COSEC I ) = <E*COS(TH( I 5 ))/(.!, +E*COS ( Tr < I > ) ) 

ECA( I)aARTAN(SINE< I ),COSE( I >,i> 
XM(I)sECA(n-E*SIN(ECA(I>) 


XMTRSXMC2) 

ETlsECA(l) 

Et2=ECA(2) 


T=SQRT( A**3/GM) 

TFA=T*XM(1) 

TFB=T*XM(2) 

IF (TFB-TFA) 20,30,30 
TFB=TFB*2 , *P I *t 
TF=TF8-TFA 


RETURN 

END 
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N 1 

N 
N 
N 

N 5 
N 6 
N 7 
N 8 
N 9 
n 10 

N 11 
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N 17 
N 18 
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N 20 
N 21 
N 22 
N 23 
N 24" 
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3200 FORTRAN ( 3 . 0 > /RTS / / 

SUBROUTINE RANGA(xT, vDT, vOMEGA,PNlT) 

DIMENSION TEMPI. <3 ) , TfcMP2C3) »TEMP3( 3) , XOMEGA (3) , XT<3) , XDT(3> 

CALL VCROSS (TEMPI, XT, XDT) 

CALL VCROSS ( TEMP2 , TEMPI , XOMEGA ) A 119 

CALL VCROSS ( TEMP3, TEMPI, TEMP2) A 120 

RRTsVMAG(XT) A 121 

SlNPMTaVD0T(TEMP3»XT)/(VMAQ(TEMP3)*RRT> A 122 

C0SPhTsVD0T(TEMP2,XT>/(.VMaQ(TeMP2)*RRT> A 123 

PHIT*ARTAN(SINPHT,C0SPHT,1> a 124 

RETURN 
END 


FORTR a N DIAGNOSTIC RESULTS FOR RaNG a 
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3200 FORTRAN ( 3 . 0 ) /RTS 



TRUE(XP.XDP,TEMPl,PTA) 
|V XP<3 > , XDP ( 3) , TEMPI( 3 ) , T 


EP = V.HAG( TEMPI ) 
RMPa VflAG ( XP ) 


EMP?(3) , TEMP3 ( 3 ) 


CALL VCROSS ( TEMP2, XDP, XP) 

CALL VCROSS C TPMP3 , TEMPI , TEMP? > 

COSPTA = v POT ( TEMPI, XP)/(EP*RMP) 
SINPtA-V0QT(T'fmP3,XP)/(VMAG<TEMP3)*RMP) 
P T A = A R T A N ( S l N P T a , C 0 S P T A , 1 ) 

R E T 0 K N 
FND 


/ 


/ 


A 549 
A 550 
A 551 
A 552 
A 553 


FORTRAN DjArNOSTjC RESULTS pOR TRU E 
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